'Si- 



Submitted to The Astrophysical Journal - revised version 

The nucleosynthesis of 26 A1 and 60 Fe in solar metallicity stars 
extending in mass from 11 to 120 M : the hydrostatic and 

explosive contributions 



Marco Limongi 1 ' 3 and Alessandro Chiem 2,3 



O 
O 
(N 

9^- ABSTRACT 

< 

We present the 26 A1 and 60 Fe yields produced by a generation of solar metal- 
licity stars ranging in mass between 11 and 12OM . We discuss the production 
sites of these 7 ray emitters and quantify the relative contributions of the various 
components. More specifically we provide the separate contribution of the wind, 



the C convective shell and the explosive Ne/C burning to the total Al yield per 



each stellar model in our grid. We provide the contributions of the He convective 
shell, the C convective shell and the explosive Ne/C burning to the 60 Fe yield as 
well. From these computations we conclude that, at variance with current beliefs, 
26 Al is mainly produced by the explosive C/Ne burning over most of the mass 
interval presently analyzed while 60 Fe is mainly produced by the C convective 
"£0 \ shell and the He convective shell. 

By means of these yields we try to reproduce two quite strong observational 
constraints related to the abundances of these nuclei in the interstellar medium, 
i.e. the number of 71. 8 photons per Lyman continuum photon, Rg x l, and the 
60 Fe/ 26 Al 7-ray line flux ratio. RqxL is found to be roughly constant along the 
galactic plane (Knodlseder 1999) (and of the order of 1.25 x 1CT 11 ), while the 
60 Fe/ 26 Al ratio has beed recently measured by both RHESSI (0.17 ± 0.05) and 
SPI(INTEGRAL) (0.11±0.03). We can quite successfully fit simultaneously both 
ratios for a quite large range of exponents of the power law initial mass function. 
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We also address the fit to 7 2 Velorum and we find that a quite large range of 
initial masses, at least from 40 to 60M Q , do eject an amount of 26 A1 (through 
the wind) compatible with the current upper limit quoted for this WR star: 
such a result removes a longstanding discrepancy between the models and the 
observational data. 

Subject headings: nuclear reactions, nucleosynthesis, abundances - stars: evolu- 
tion - stars: supernovae 



1. Introduction 

The quest for the main 26 A1 source (s) started as soon as the first observational evidence 
of the presence of live 26 A1 in the inner Galaxy was demonstrated by the HEAO 3 experiment 
(Mahoney et al. 1982, 1984) and it still continues at the present days. Its presence in the 
interstellar medium was not a surprise because it was theoretically predicted since the middle 
of the seventies that explosive Ne/C burning in core collapse supernovae may produce such 
a nucleus (Arnett 1977; Ramaty & Lingenfelter 1977; Truran & Cameron 1978; Arnett & 
Wefel 1978; Woosley & Weaver 1980). Later on Dearborn & Blake (1985) pointed out that 
26 Al may by injected into the interstellar medium also by Wolf-Rayet stars through stellar 
wind and since then a number of works have been devoted to the prediction of the amount 
26 Al ejected by these stars (e.g. Prantzos & Casse 1986; Walter & Maeder 1989; Meynet et 
al. 1997; Vuissoz et al. 2004; Palacios et al. 2005). Arnould et al. (1980) and Clayton (1984), 
on the basis of the theoretical models available at that time, proposed also the novae as main 
26 Al producers. All these candidates, with the further addition of the AGB stars (Mowlavi 
& Meynet 2000), could in principle account for a significant fraction of the observed flux due 
to the uncertainties in both the frequencies of the various candidates and the yields itself. 
Leising & Clayton (1985) were the first to point out that the lack of information about the 
distribution of the 26 A1 in the Galaxy severely limits the possibility to identify the main 26 A1 
donor (s). A clear snapshot of the situation just before the launch of the Compton Gamma 
Ray Observatory (CGRO) is provided by the Prantzos (1991) review paper. 

The launch of the CGRO in 1991, and its continuous operation up to 2001, allowed 
the first mapping of the 1809 KeV line all over the Galaxy and put firm constraints on 
the possible main sources of 26 A1. By the way, the latest all-sky map has been published 
by Pliischke et al. (2001). This observational campaign showed beyond any doubt that 
the 26 A1 distribution was strongly confined towards the galactic disk and quite clumpy and 
asymmetric, strongly pushing towards a massive star parent for the bulk of the 26 A1. We 
refer the reader to the reviews by Prantzos & Diehl (1996) and Diehl & Timmes (1998) for 
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an excellent clear picture of the scenarios that were under debate as soon as the first results 
from COMPTEL aboard the CGRO became available. A further milestone progress in the 
identification of the main 26 A1 producers was obtained when the 26 A1 map was compared 
to a series of all sky maps obtained at various wavelengths. Knodlseder et al. (1999) and 
Knodlseder (1999) found that the 53 GHz microwave free- free emission all-sky map was the 
only one to almost perfectly overlap the 26 A1 all-sky map. The 53 GHz free-free emission 
line maps the regions of ionized gas and such a sustained ionization may be powered only 
by the Lyman continuum photons provided by the very massive stars and disappears on a 
timescale of 1 Myr after the switch off of the Lyman continuum "lamp" . Hence the existence 
of such a good overlap between the two all-sky maps demonstrates beyond any reasonable 
doubt that the bulk of the 26 A1 is produced in the same regions where massive stars are 
present. Such a finding did not close at all the quest of the 26 A1 progenitors but it simply 
shifted such a quest towards which subclass of massive stars do mainly contribute to the 26 A1 
budget. The two classical competing candidates were identified as the Type II supernovae 
and the Wolf Rayet (WR) stars (through the stellar wind). A further question concerned 
the possibility, sometimes suggested, that a consistent fraction of the 1.8 MeV 7-flux (71. s) 
were due to close localized sources and not to Galaxy wide sources. Such a doubt was 
mainly due to the fact that the massive star census in the Galaxy is well known only out 
to distances of a few kiloparsecs and many regions of the Galaxy are obscured for direct 
measurements. This problem was recently definitely solved by SPI(INTEGRAL) (Diehl et 
al. 2006) that measured for the first time with great accuracy the doppler shift of the 26 A1 
line along the inner Galaxy clearly demonstrating that the 26 A1 corotates with the Galaxy, 
hence supporting a Galaxy-wide origin. The total amount of 26 A1 present in the disk of our 
Galaxy has been estimated by assuming a variety of distributions of the 26 A1 sources and it 
ranges between 1 and 3 M Q . 

Together with the search for a signal from the decay of 26 A1, also a signal from an- 
other gamma ray emitter, namely 60 Fe, has been searched for thoroughly. The production 
of this nucleus has been historically attributed to the Ne explosive burning in core collapse 
supernovae (Clayton 1982), hence it has always been considered as the litmus paper to 
discriminate among the various candidates as major 26 A1 contributors. Unfortunately, all 
the experiments up to 2002 could only provide upper limits on the abundance of this nu- 
cleus in the interstellar medium. With the launch of the Reuven Ramaty High Energy Solar 
Spectroscopic Imager (RHESSI) and of the International Gamma-Ray Astrophysics Labo- 
ratory (INTEGRAL) the situation changed drastically since both RHESSI and SPI aboard 
INTEGRAL were able to measure an 60 Fe signal from the inner radiant of our Galaxy. The 
latest simultaneous measurements of 26 A1 and 60 Fe performed by RHESSI and SPI at the 
time of writing give a 60 Fe/ 26 Al 7-ray line flux ratio equal to (0.17 ± 0.05) (RHESSI, Smith 
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2005) and (0.11 ±0.03) (SPI, Harris et al. 2005), respectively. 

The theoretical models that constitute the base of the various interpretations of these 
data may be divided in two broad categories: the first one includes stars less massive than, 
say, 35-40M Q and a second one that includes the more massive stars. This limiting mass 
roughly marks the minimum mass that becomes a WR star. The evolution of the stars in 
the first group has been extensively explored through all the hydrostatic burnings up to 
the core collapse together with the following explosive burnings. Yields for both 26 A1 and 
60 Fe have been published, for this mass range, over the years by several groups (Woosley & 
Weaver 1995; Thielemann et al. 1996; Rauscher et al. 2002; Limongi & Chieffi 2003; Chieffi 
& Limongi 2004). The evolution of the stars in the upper mass range, which spend part 
of their lifetime as WR stars, have never been followed all along their full evolution (and 
final explosion) but only through the H and He burning (plus, quite recently, the central C 
burning) phases; obviously also the explosive burnings of these very massive stars have never 
been addressed. As a consequence of this very biased knowledge of the evolution of these very 
massive stars, only the fraction of 26 A1 yield ejected by the stellar wind has been predicted 
up to now (e.g. Prantzos et al. 1986a,b; Walter & Maeder 1989; Meynet et al. 1997; Palacios 
et al. 2005), and unfortunately quite often this partial prediction has been interpreted as 
the total 26 A1 yield produced by these very massive stars. In other words, it has become a 
deep-seated notion that the stars that do not become WR contribute to both the 26 A1 and 
60 Fe yields while those that become WR do contribute only to the synthesis of the 26 A1. 
Within this accepted framework, an "observed" 60 Fe/ 26 Al flux ratio of 0.11-0.17 would favor 
the masses that do not become WR as the main producers of both 26 A1 and 60 Fe (Timmes 
et al. 1995). Indeed, the WW95 models do predict a 60 Fe/ 26 Al flux ratio of 0.18 (averaged 
over a power law IMF - dn/dm = M^ 235 ^ - ranging between 11 and 40M Q ). However, more 
recent models (Rauscher et al. 2002; Limongi & Chieffi 2003; Chieffi & Limongi 2004) do 
predict much larger 60 Fe/ 26 Al flux ratios (0.5- 1.0), leaving therefore plenty of room for an 
important contribution of the WR stars to the synthesis of the 26 Al. For sake of completeness, 
it must be mentioned that Woosley et al. (1995) published the presupernova evolution and 
explosive yields of massive pure He cores that should be associated to WR stars. However, 
unfortunately, these bare cores can hardly be reconciled with a population of (single) WR 
stars because all of them correspond to models that lose the H-rich mantle before the central 
He ignition while consistent models of massive stars with mass loss (see below) show that 
this happens only in stars more massive than 80 M Q . 

In this paper we completely reanalyze and rediscuss the interpretation of the observed 
fluxes of both 26 A1 and 60 Fe in the light of a new large grid of massive star models of solar 
metallicity with mass loss. These homogeneous sets of models extend in mass between 11 and 
12OM and cover the full hydrostatic evolution from the Main Sequence up to the core col- 
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lapse together with the following explosive burnings. These models will be published shortly. 
Given the importance of the problems connected with the 7 ray emitters 26 A1 and 60 Fe, we 
think that it is worth dedicating a full paper to them. We are aware that these models are 
just solar metallicity models and that it would be important to have supersolar models too, 
but we think that the availability of a full set of models from 11 to 120M Q covering all the 
hydrostatic and explosive burnings is already a big step forward in the understanding of 
these nuclei. We also address the fit to 7 2 Velorum, the closest WR star. 

The paper is organized as follows: first, the new set of models is briefly presented, then 
two distinct sections discuss separately the theoretical yields of these two nuclei as a function 
of the initial mass. Their dependence on the mass loss rate adopted in the WNE+WCO 
phase is addressed in Section 5. A discussion of the scenarios that comes out from the 
adoption of these new yields and a final brief conclusion follow. 



2. The stellar models 

The 26 A1 and 60 Fe yields used in the present paper have been extracted from a new gen- 
eration of stellar models that will be available soon. These new models have been computed 
with the latest release (5.050218) of the FRANEC (Frascati RAphson Newton Evolutionary 
Code). By the way, from now on the release number will follow the following convention: 
the integer part of the number refers to the major release (the fifth at present) while its 
decimal part refers to the date in which it has been released, in the form YYMMDD. The 
main features of this new release, with respect to that described in Limongi & Chiefn (2003), 
are the following. First of all the convective mixing and the nuclear burning are coupled 
together and solved simultaneously. More specifically, the set of equations describing the 
chemical evolution of the matter are now: 



m = f&Yf\ + _d_ 
dt V dt / nuc dm 



am 



where the diffusion coefficient D is given by D = l/3vj, the convective velocity v c 
is computed by means of the mixing length theory and I is the mixing lenght parameter. 
This set of equations is linearized and solved by means of a Raphson Newton technique. 
Let us also remind that, as in all our computations, the nuclear energy generation rate is 
always computed with the same, full, network adopted to describe the chemical evolution 
of the matter, without any kind of approximation. The network adopted in the present 
computations is the same as in Limongi & Chieffi (2003). 0.2 H p of overshooting has been 
assumed at the top of the convective core in central H burning. The nuclear cross sections 
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have been updated with respect to those adopted in Limongi & Chieffi (2003) whenever 
possible. Table 1 shows the full reference matrix of all the processes taken into account 
in the network, together with its proper legend. Note that for the weak interaction rates, 
f3 + and j3~ mean the sum of both the electron capture and the f3 + decay and the positron 
capture and the f3~ decay, respectively. The dependence on the temperature of the weak 
interaction rates has been considered following Oda et al. (1994), Fuller, Fowler & Newman 
(1982), Langanke & Martinez-Pinedo (2000) and Takahashi & Yokoi (1987), see Table 1. The 
26 Al requires a special treatment since the ground (Alg6 entry in Table 1) 26 Al g (J 71 = 5 + , 
r = 1.03x 10 6 yr) and the isomeric states (Alm6 entry in Table 1) 26 Al m (.F = + , r = 9.15 s) 
do not come into statistical equilibrium for temperatures T < 10 9 K (Ward & Fowler 1980; 
Gupta & Meyer 2001). For this reason we treat the ground and isomeric states of 26 A1 as 
separate species for T < 10 9 K while we assume the two states to be in statistical equilibrium 
(and hence consider just one isotope) above this temperature. By the way, the A126 entry 
in Table 1 refers to the total amount of 26 A1. 

Mass loss has been included following the prescriptions of Vink et al. (2000, 2001) for the 
blue supergiant phase (T c g- > 12000 K), de Jager et al. (1988) for the red supergiant phase 
(T cff < 12000 K) and Nugis & Lamers (2000) for the Wolf-Rayet phase. Additional models 
have been computed with the Langer (1989) mass loss rate in the WNE+WCO phases. 

The explosion of the mantle of the star is started by imparting istantaneously an initial 
velocity Vq to a mass coordinate of ~ 1 M Q of the presupernova model, i.e., well within 
the iron core. The propagation of the shock wave is computed by means of an explosive 
simulation code, developed by us, that solves the fully compressible reactive hydrodynamic 
equations using the Piecewise Parabolic Method (PPM) of (Colella & Woodward 1984) in the 
lagrangean form. Vq is properly tuned in order to eject all the mass above the Fe core. The 
chemical evolution of the matter is computed by coupling the same nuclear network adopted 
in the hydrostatic calculations to the system of hydrodynamic equations. The nuclear energy 
generation is not taken into account in these simulations since we assume that it is always 
negligible with respect to both the kinetic and the internal energies. By taking advantage of 
the fact that the final yields mainly depend on the mass cut location (the mass coordinate 
that separates the ejecta from the compact remnant) and only mildly on the properties of 
the parametrized explosion, yields corresponding to different amounts of 56 Ni ejected can be 
easily obtained by fixing the mass cut by hand a posteriori (Limongi & Chieffi 2003). Here, 
for each star, we chose the mass cut in order to obtain 0.1 M of 56 Ni ejected. By the way, 
let us note that since both 26 A1 and 60 Fe are synthesized in quite external regions, their final 
yields do not depend on the particular choice of the mass cut since we assume in any case 
that at least a minimum amount of 56 Ni is ejected. 
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All these solar models has been computed by adopting an initial He abundance equal 
to 0.285 (by mass fraction), a global metallicity Z=0.02 (by mass fraction) - the relative 
abundances among the various nuclear species are taken from the Anders & Grevesse (1989) 
solar distribution. They range in mass between 11M and 120M Q , covering therefore the 
full range of masses that are expected to give rise to the Type II/Ib/Ic Supernovae as well 
as to the WR populations. These models and the corresponding full set of yields will be 
published, and discussed, shortly; for the moment we anticipate a few basic properties in 
Table 2. The following data are reported: the initial and final mass in solar masses (cols. 1 
and 2), the central H and He burning lifetimes in Myr (cols. 3 and 4), the lifetime after the 
central He exhaustion in yr (col. 5), the maximum sizes of the He and CO core masses in 
solar masses (cols. 6 and 7), the bottom and top mass locations of the last C convective shell 
in solar masses (col. 8), the central C mass fraction at the central He exhaustion (col. 9), the 
total amount of time spent as O star (col. 10) and WR star (col. 11) in yr and eventually 
the lifetimes as WNL (col. 12), WNE (col. 13) and WCO (col. 14) all in yr. We assume that 
the star enters the WR phase when Log(T c g) > 4 and H sur f < 0.4 and we adopt the following 
usual definitions for the various WR phases: WNL (10~ 5 < H surf < 0.4), WNE (H surf < 10~ 5 
and (C/N) surf < 0.1), WNC (0.1 < (C/N) surf < 10) and WCO ((C/N) surf > 10). The main 
data concerning the 26 A1 and 60 Fe yields are summarized in Table 3. The first column refers 
to the initial mass in solar masses. Columns 2 to 5 refer to the 26 A1 yield (all in solar masses): 
total (col. 2), the wind contribution (col. 3), the C&Ne shell contribution (col. 4) and the 
explosive component (col. 5). Columns 6 to 10 refer to the 60 Fe yield (all in solar masses): 
total (col. 6), the contributions of the radiative and convective He burning shell (cols. 7 to 
9), the convective C shell contribution (col. 10) and the explosive component (col. 11). 



3. 26 A1 

26 A1 is an unstable nucleus (its terrestrial half life is r\ii — 7.17 x 10 5 y) produced almost 
exclusively by proton capture on 25 Mg. Its synthesis occurs, in massive stars, essentially in 
three different specific environments, i.e., the core H burning, the C and Ne convective shells 
and the explosive Ne burning. Let us discuss each of these environments separately. 

All over the mass range presently analyzed, 26 A1 is produced in Main Sequence by the 
25 Mg(p,7) 26 Al reaction (hereinafter MG25PG) and mainly destroyed by the f3 + decay into 
26 Mg since the competing destruction process, i.e. the 26 Al(p,7) 27 Si process, becomes efficient 
for T > 50 x 10 6 K and only our most massive model barely reaches such a temperature 
towards the end of the central H burning. 

The conversion of the initial abundance of 25 Mg into 26 A1 starts as soon as a star settles 



on the Main Sequence (since T centra i > 30 x 10 6 K) and comes to completion on a timescale of 
few million years, i.e. well within the central H burning lifetime. Figure 1 shows the temporal 
evolution of the 26 A1 mass fraction in the H convective core: a maximum is reached shortly 
after the beginning of the central H burning, followed by a progressive decrease mainly due 
to the f3 + decay. Since 26 A1 continues to be produced also after the maximum, its decline is 
shallower than that of a pure decay (shown in Figure 1 as the cyan dot-dashed line). By the 
way, it is such a prolonged production that preserves some 26 A1 in the H convective core up 
to the end of the central H burning. In fact, if the slope of the decline would follow a pure 
decay after the maximum, essentially all the 26 A1 would decay by the time the star exhausts 
the H in the center. The solid, long dashed, dotted and short dashed lines refer to the 15, 
30, 60 and 120 M© , respectively. 

Since the efficiency of the MG25PG reaction rate scales directly with the initial mass, 
the larger the initial mass the higher the maximum 26 A1 mass fraction abundance, the faster 
the 25 Mg consumption and hence the faster the 26 A1 decline. Hence, if the efficiency of the 
MG25PG were the only parameter to control the scaling of the 26 A1 mass fraction abundance 
at the central H exhaustion with the initial mass of the star, one should obtain an inverse 
scaling of the central 26 A1 mass fraction at the central H exhaustion with the mass. Vice 
versa the models show a strong direct scaling. Such a result is the obvious consequence of 
the strong shortening of the central H burning lifetime with the initial mass that freezes the 
central 26 A1 abundance progressively at earlier times (and hence at higher abundances). Of 
course, a direct scaling between the 26 A1 mass fraction and the initial mass of a star does not 
imply, by itself, a direct scaling between the total amount of 26 A1 (in solar masses) produced 
by the H burning and the initial mass. In fact, it is the size of the convective core that really 
determines how much 26 A1 is produced in total, and since it scales directly with the initial 
mass, the direct scaling between 26 A1 produced and initial mass is secured. 

At central H exhaustion, the 26 A1 is located both in the He core and in the region of 
variable H left behind by the receding convective core. Since the He burning easily and 
quickly destroys the 26 A1 (via the (n, a) and (n, p) reactions - the neutrons being released 
by the l3 C(a, n) 16 process) the amount of 26 A1 synthesized by central H burning and really 
preserved up to the explosion is just the one located in the H rich layers plus the one locked 
in the fraction of the He core that will not be affected by the He burning. If the dredge up 
and the mass loss were not effective, this 26 A1 would be ejected all together by the explosion, 
while the occurrence of these two phenomena may anticipate such an ejection at earlier times 
and also increase the amount of 26 A1 ejected into the interstellar medium. In our models 
three mass intervals may be identified: 1) stars between 11 M and 30 M undergo a dredge 
up episode that does not enter into the He core and the mass loss is weak enough that all 
these stars end their life with an H-rich envelope; since the He convective shell extends almost 
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up to the base of the H burning shell, only the tiny amount of 26 A1 present in the region of 
variable H left by the receding H convective core and engulfed in the convective envelope is 
ejected in the interstellar medium. 2) stars between 35M Q and 4OM experience a dredge up 
episode that does not enter the He core as well, but the mass loss is strong enough that also 
a fraction of the He core is ejected outward; since the mass loss erodes the He core before 
the formation of the He convective shell, the 26 A1 present in these layers is preserved by the 
destruction and goes to increase the amount of 26 A1 ejected in the interstellar medium. 3) 
stars more massive than 4OM do not show any dredge up episode and the mass loss is so 
strong that a substantial fraction of the He core is ejected through the stellar wind; as a 
consequence a larger amount of 26 A1 present in the He core is preserved from the destruction 
and ejected. As for the first two mass intervals, the wind starts expelling 26 A1 as soon as 
the star becomes a red (super)giant. Vice versa, stars in the upper mass interval begin to 
eject 26 A1 only when the total mass of the star reduces enough that layers processed by the 
H burning are exposed to the surface (i.e. when the star becomes a WR). Let us eventually 
note that only a very tiny fraction of the original 25 Mg is ejected in the interstellar medium as 
26 Al(mainly because of the substantial 26 A1 decay during the core H burning) and also that 
this 26 A1 is of semi- secondary origin since only the 25 Mg depends on the initial metallicity Z 
and obviously not the protons (the true secondary nuclei are the ones formed by reactions in 
which both the target and the projectile depend on the initial metallicity, for example, as in 
the s process production). Figure 2 shows, as green diamonds, the amount of 26 A1 expelled 
through the stellar wind as a function of the initial mass; the observed trend is the obvious 
consequence of what has been discussed above, i.e., it is the result of the direct scaling of 
(1) the size of the H convective core, (2) the 26 A1 mass fraction left by the central H burning 
and (3) the mass loss efficiency with the initial mass. 

Once the star moves beyond the central H burning, the "habitat" becomes very hostile 
to the survival of 26 A1. First of all, how already mentioned above, the He burning simply 
destroys 26 A1 through neutron captures. Second, the half life of 26 A1 has a tremendous 
dependence on the temperature: it reduces to ~ 0.19 yr already at Log(T) = 8.4, to ~ 13 
h at Log(T) = 8.6, and to 13 4- 2 m in the range 9 < Log(T) < 9.6. Third, efficient 
photodisintegration quickly destroy 26 A1 when the temperature raises above, say, Log(T) = 
9.3. Such a temperature limit obviously constrains a possible 26 A1 production to the C 
and/or Ne burning shells. 

As for the C shell, its typical temperature (Log(T) < 9.08) does not allow in general a 
substantial production of 26 A1. However, after the central Si burning, the strong contraction 
(and heating) of the inner core that precedes the final gravitational collapse, induces a strong 
temperature increase of the C burning shell (up to, roughly, Log(T) ~ 9.255). If, at this 
stage, such a burning occurs in an efficient convective shell, a substantial amount of 26 A1 is 
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produced. Let us discuss such a production in a more detailed way. The processes that pro- 
duce and destroy 26 A1 in these conditions are once again the MG25PG and the beta decay. 
The MG25PG reaction rate depends on the MG25PG cross section but also on the abun- 
dances of 25 Mg and protons. The 25 Mg that now enters in the 26 A1 production comes directly 
from the initial CNO abundance, via the sequence (CNO)j n i — » 14 N— » 22 Ne— > 25 Mg while the 
protons are mainly produced by the 12 C( 12 C,p) 23 Na and the 23 Na(o;,p) 26 Mg processes. By 
the way, note that the main proton poison is not the 25 Mg but, instead, the 23 Na through 
the (p,a) channel. The high temperature influences the 26 A1 production either because it 
raises the MG25PG cross section but also because it strongly increases the proton density by 
increasing the cross sections of both the processes that provide the bulk of protons. On the 
other hand, the presence of a convective environment is also crucial because it continuously 
brings fresh 12 C and 25 Mg where the 26 A1 production occurs and simultaneously brings the 
freshly made 26 A1 to much lower temperatures where its lifetime increases enormously. 

At variance with the C burning, 26 A1 is always produced by the Ne burning (both 
in the center and in the shell). A substantial amount of 25 Mg is left unburned by the 
C burning and the protons needed to activate the 25 Mg(p,7) 26 Al reaction mainly come 
from the 23 Na(a,p) 26 Mg reaction ( 23 Na being one of the products of C burning), the oc's 
being provided by the 20 Ne(7, a) 16 photodisintegration. Again, the most effective protons 
poison is the 23 Na(p, a) 20 Ne. In a radiative environment the 26 A1 "equilibrium" abundance 
depends on the local balance between production and destruction so that, as soon as the 
20 Ne photodisintegration weakens for lack of fuel, the abundance of the a's drops as well 
and the proton production channels dry up: no more 26 A1 may be produced and its short 
lifetime at these high temperatures leads quickly to its total destruction. Once again only the 
presence of an efficient convective shell could act as a preservation buffer. Since most of the 
stars do not reach the onset of the explosion with an efficient Ne convective shell, in general 
only the small equilibrium abundance of 26 A1 located within the radiative Ne burning shell 
is present at the beginning of the core collapse. There are two exceptions: the 14 and the 15 
M Q models. Stars within this mass interval are characterized by the lack of a C convective 
shell, at the end of the central Si burning, and by the formation of a Ne convective shell 
very close to the region where it was previously efficient a carbon convective shell. During 
the last strong contraction of the core that precedes the final collapse, such a Ne convective 
shell penetrates the C rich region with the consequence of producing a huge amount of 26 A1. 

The 26 A1 synthesized in the C (or Ne/C) convective shell is located close enough to 
the iron core that is partially destroyed by the passage of the shock wave. In general, the 
larger the mass the larger the amount of 26 A1 that survives the explosion (Table 3). The blue 
triangles in Figure 2 show the 26 A1 yield produced in this convective shell that survived to the 
explosion as a function of the initial mass. The general increase of this yield with the mass 
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depends on the fact that both the mass size and the peak temperature of the C (or Ne/C) 
convective shell scale directly with the mass. The peak around the 15M model, as already 
discussed above, is due to the presence of an extended Ne convective shell that engulfs part 
of the C rich region. The minimum yield around the 3OM is the natural consequence of 
the lack of an active C (and Ne) convective shell at the moment of the explosion. The 26 A1 
produced by the C and/or Ne shell is of semi-secondary origin (as the one produced by the 
H burning) because the 25 Mg is secondary while protons are primary. 

26 Al is also produced during the explosion at a typical temperature of the order of ~ 2.3 
billion degrees. Such a condition occurs within the C convective shell and the main process 
that controls its production is once again the MG25PG process while its destruction is now 
controlled, roughly parithetically, by the (n, p) and (n, a) processes. The 25 Mg now comes 
mainly from the (n, 7) capture on 24 Mg, this isotope being a primary outcome of the C and 
Ne burning. The neutron density that enters in both the production and destruction of 26 A1 
is determined by the competition among several processes: the main neutron producers are 
(a, n) captures on 26 Mg, 25 Mg, 21 Ne, 29 Si plus the (p, n) capture on 28 A1 while the main 
neutron poisons are the (n, 7) captures on 24 Mg, 16 and 20 Ne. The proton density that 
enters the MG25PG rate is mainly determined by the competition between the production 
that occurs via the (a, p) captures on 20 Ne, 24 Mg, 27 A1 and 23 Na and the destruction that 
occurs via the (p, 7) capture on 26 Mg, 20 Ne, 24 Mg, 27 A1, 25 Mg, 30 Si plus the (p, n) reaction on 
28 Al. Figure 2 shows the contribution of the explosive burning to the synthesis of 26 A1 as filled 
red squares. To understand the trend of the explosive yield with the mass it is necessary to 
remind that the final Mass Radius (MR) relation at the moment of the core collapse plays a 
fundamental role for the synthesis of many nuclei, because it mainly determines the amount 
of matter exposed to the various explosive burnings (Chieffi et al. 2000; Chieffi & Limongi 
2002). Such a relation strongly depends on the behavior of the C, Ne and O shell burnings 
all along the evolution of a star. The general rule is that the stronger is a shell burning, 
the slower is the contraction and hence the shallower is the final MR relation; it is also 
worth reminding that a stronger shell burning usually implies a wider (in mass) convective 
shell. Figure 3 shows the final MR relation for a subset of stars in the present grid. By the 
way, since the blast wave is essentially radiation dominated, the spatial location where the 
peak temperature drops to, say, 2.3 billion degrees occurs around 0.01-0.02 R , the exact 
value depending (quite mildly) on the explosion energy. The Figure clearly shows that there 
is a non monotonic dependence on the initial mass: though a general trend exists, in the 
sense that the larger the mass the steeper this relation, two inversions exist: the first one 
between 17M Q and 2OM and the second one between 3OM and 35M . The first inversion 
marks the first mass (2OM ) forming a very extended last C convective shell that slows 
down the gravitational contraction of the core and hence halts the direct scaling between the 
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compactness of the inner layers and the initial mass. The second inversion, occurring around 
the 35M model, marks the minimum mass that forms a single well settled C convective 
shell that lasts up to the explosion. Also in this case the formation of a single very active 
C convective shell counterbalances the increase in the gravitational contraction induced by 
the increase of the total mass. The explosive yields shown in Figure 2 follow almost exactly 
the general trend shown by the MR relation: there is a general direct scaling between the 
explosive 26 A1 and the initial mass, but with two minima corresponding to the two masses 
quoted above. 

The black dots in Figure 2 show the total yield of 26 A1 as a function of the initial 
mass. It is quite evident that for the largest majority of the stars in the present mass range 
11 < M/M < 120 the explosive component is the major contributor to the total yield. Only 
around 15M the 26 A1 produced in the advanced hydrostatic burnings dominates the final 
yield. The wind component to the total is always negligible up to a mass of, say, 6OM . 
Above this mass, it contributes to roughly 30% of the global outcome. Note, however, that 
though the explosive burning is the overall major producer of 26 A1 within the full initial mass 
interval presently analyzed, the 26 A1 produced by both the H burning and the C/Ne shell 
cannot be ignored in the computation of the final budget of 26 A1 produced by a generation 
of massive stars. 

There are many quite uncertain basic stellar parameters that may significantly alter the 
final 26 A1 yields shown in Figure 2. We have explored a few of them: 

1) the adopted initial 25 Mg abundance (usually scaled solar) constitutes the main buffer 
from which the 26 A1 is produced in MS and a change of its initial (assumed) abundance 
reflects linearly on the final amount of 26 A1 ejected by the wind. In other words, a doubling 
of the initial 25 Mg leads to a doubling of the 26 A1 ejected through the wind. 

2) the size of the H convective core certainly plays a pivotal role since it affects both 
directly and indirectly the final 26 A1 yield. The direct influence occurs in MS because it 
contributes to determine its abundance, by mass fraction, as well as its distribution within 
the star at the end of the core H burning; the indirect influence occurs through its major 
role in determining the He core mass, which is the main parameter driving all the further 
evolution of a star. In order to study the dependence of the 26 A1 yield on the size of the H 
convective core we recomputed the full presupernova evolution plus the explosive burning of 
a 60 and a 12OM by adopting 0.5 Hp of overshooting (instead of the basic 0.2 Hp) in core H 
burning. The main results are shown in the last two rows of Table 2 and Table 3. The first 
thing worth noting is that the increase of the size of the H convective core has two opposite 
effects on the two masses: the total 26 A1 yield increases by 18% in the 6OM while it remains 
roughly constant in the 12OM . A closer look to the various components shows that the 
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wind component reduces by 16% in the 60M Q while it increases by 24% in the 120M Q . 
The sum of the other two components shows, on the other hand, the opposite trend, i.e. it 
increases in the smaller mass and reduces in the larger one. Such results are the consequence 
of a complex interplay among several phenomena. First of all an increase of the convective 
core has two opposite effects: it increases the region where the initial 25 Mg is converted in 
26 Al but it also leads to a lower 26 A1 abundance by mass fraction at the core H exhaustion 
(because the H-burning lifetime increases). The final balance between these two opposite 
phenomena is not predictable a priori, but it happens that the first effect prevails on the 
second one so that at last more 26 A1 (in solar masses) is produced overall. However, such an 
occurrence is not sufficient to determine how the 26 A1 wind component depends on the size of 
the convective core because another crucial role is played by mass loss: a change of the size of 
the convective core strongly affects the path followed by a star in the HR diagram and hence 
the total amount of mass it loses. In the 6OM case, the standard and test models lose similar 
amounts of mass in MS; however, while the standard model spends a consistent fraction of 
the central He burning phase as a Red Supergiant where it loses a very large amount of mass, 
the model with the larger convective core moves much earlier to the blue where the average 
mass loss rates are significantly lower. The final result is that the test model loses less mass, 
reaching the final collapse with a total mass of 19.93M , i.e. 2.8M more than the standard 
case. This explains why the test 6OM ejects less 26 A1 through the wind than its standard 
counterpart. In the 12OM case the opposite happens: the test model loses much more mass 
than the standard model: in fact its final mass is 26.83M , i.e. 4M less than the standard 
model. In this specific case the luminosity increase in central H burning (due to the increase 
of the convective core) allows the model to enter (for a consistent amount of time) a region of 
the HR diagram where mass loss is very efficient (see Vink et al. 2000, 2001). This explains 
why the amount of 26 A1 present in the wind is larger in the test than in the standard model. 
As for the 26 A1 produced by the advanced (C/Ne) burnings plus the explosion, we found 
that these components scale directly with the size of the convective core in the 6OM and 
inversely in the 12OM . The reason is simply that the amount of 26 A1 produced both in the 
advanced and in the Ne explosive burnings scale directly with the size of the He core mass, 
and Table 2 shows that this parameter increases with the size of the convective core in the 
6OM , while it decreases in the 12OM case: this last (apparently) unexpected result is due 
to the strong increase of the mass loss rate in the 120 M stellar model computed with 0.5 
Hp of overshooting. 

3) We have also explored the influence of the MG25PG cross section on the 26 A1 produced 
by the H burning and expelled through the wind. Though the uncertainty quoted for this 
process (e.g. Iliadis et al. 2001) is of the order of a factor of 2 in the temperature range 
0.02 < T9 < 0.2, we preferred to explore the general theoretical dependence of the 26 A1 
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(ejected by the wind) on the MG25PG cross section in a range of values wider than (presently) 
expected. Hence we have performed a few tests on both the 60M Q and the 12OM by 
multiplying and dividing this cross section by a factor of 3 and 10. Table 4 shows the results 
of these tests in the first four rows: the first column shows the cross section adopted in the 
test while the forth and fifth columns show the amount of 26 A1 (in solar masses) present 
in the wind for the 6OM and the 120M Q star models, respectively. These tests show that 

a) the amount of 26 A1 ejected by the wind does not scale monotonically with the efficiency 
of this cross section but, instead, that it firstly increases as the cross section decreases and 
then it begins to decrease when the cross section drops below 1/3 of the reference one, and 

b) the amount of 26 A1 ejected by the wind depends only mildly on a systematic variation 
of the cross section, since it varies by less than a factor of two over a huge variation of 
the MG25PG cross section. To understand these results it is necessary to analyze in more 
details the synthesis of 26 A1 in H burning. In the following we will discuss the behavior of the 
60M Q but the same analysis does hold over the full mass range presently analyzed. Figure 
4 shows the variation of the abundance of the 26 A1 in the convective core as a function of 
time for the five values of the cross section: the red sparse dotted line refers to the case in 
which the enhancement factor (EF) of the standard cross section is 10, the blue long dashed 
one to an EF=3, the black solid line to the standard case, the green short dashed line to an 
EF=l/3 while the magenta short dotted line refers to an EF=0.1. The test with the highest 
enhancement factor simulates the extreme case in which the 25 Mg is completely converted 
in 26 A1 just at the beginning of the H burning. In this case the maximum 26 A1 abundance is 
the highest but also the following decline is the steepest (it practically follows a pure decay) 
because no more 25 Mg is left to feed the 26 A1 abundance at late times. As the enhancement 
factor reduces, the maximum equilibrium abundance of the 26 A1 drops accordingly while the 
final 26 A1 abundance progressively increases due to the larger availability of 25 Mg at late 
time. This trend, however, naturally comes to an end when the cross section becomes so 
slow that the timescale over which 25 Mg converts in 26 A1 becomes longer than the H burning 
timescale. At this point the final abundance of 26 A1 is mainly controlled by the total amount 
of 25 Mg destroyed: the lower the rate, the smaller the amount of 25 Mg converted in 26 A1 
and hence the smaller its final abundance. A cross section 10 times lower than our standard 
case is low enough to enter such a regime. By the way, it goes without saying that the 
physical evolution of the star does not depend at all on the rate of this process, so that all 
the timescales, the amount of mass lost, the penetration of the dredge up and so on remain 
rigorously identical in all these tests. 

4) Since most of the 26 A1 ejected in the interstellar medium is anyway synthesized by the 
passage of the shock wave, we have also explored its sensitivity to the rates of the processes 
involved in its production in Ne explosive burning. Of course the number of factors that may 
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influence the explosive production of 26 A1 is very large and includes also all the factors that 
may affect the final mass-radius relation (as, e.g., the 12 C(o;,7) 16 0); such a deep analysis 
goes well beyond the purposes of this paper. We have explored the influence of a few rates, 
namely the two main processes involved in the synthesis of 26 A1, i.e. the 24 Mg(n, 7) 25 Mg 
and the 25 Mg(p,7) 26 Al, as well as the two processes that control its destruction, i.e. the 
26 Al(n, p) 26 Mg and the 26 Al(n, a) 23 Na. In particular we have recomputed three explosions 
for the two masses 25 and 6OM , the first one doubling the cross sections of both the neutron 
captures on 26 A1, a second one doubling the n capture on 24 Mg and a last one by doubling the 
MG25PG. The results are shown in rows 8 to 11 in Table 4. Since the efficiency of any process 
depends on its rate (i.e. the cross section times the abundances of the involved nuclei) rather 
than just on its cross section, the tests shown in Table 4 may also be reinterpreted as due to 
a change by a factor of two of any of the two nuclear species involved in the given process. 
This means that the tests concerning in particular the influence of the 24 Mg(n, 7) 25 Mg, the 
26 Al(n, p) 26 Mg and the 26 Al(n, a) 23 Na on the explosive synthesys of the 26 A1, may also be 
seen as tests in which it is the neutron density to have been increased by a factor of two. 
Since a doubling of the 24 Mg(n, 7) 25 Mg rate increase the explosive 26 A1 by roughly 60% 
while the simultaneous doubling of both the 26 Al(n, p) 26 Mg and the 26 Al(n, a) 23 Na rates 
reduces the 26 A1 by roughly 40%, it is clear that the two effects almost cancel out, so that 
we feel confident to conclude that a changing of the neutron density plays a minor role 
on the explosive synthesys of the 26 A1. Let us stress here, as a general warning, that the 
neutron and proton equilibrium abundances that enter in the various processes depend on 
the balance among a variety of similarly efficient processes, so that all of them should be 
carefully checked in order to understand if these equilibrium abundances are reliable or not. 

Let us eventually explicitly mention other important sources of uncertainty that could 
affect the 26 A1 yield: 

a) the C abundance left by the He burning provides the main fuel that feeds both the C 
and the Ne burning (whose abundance derives directly from the C one) and hence it regulates 
the birth and the development of all the convective shells and hence, in turn, the final MR 
relation. Unfortunately, as it is well known, this value is not firmly established because it 
depends on the still very uncertain cross section of the 12 C(o;,7) 16 process together with 
the adopted mixing scheme in central He burning. 

b) we have shown how critical the presence of a convective shell and the efficiency of 
mixing is to preserve the 26 A1 produced in the C (Ne/C) burning shell; any change in their 
mass size and/or in the efficiency of mixing could significantly alter their contribution to the 
total yield (see also Weaver & Woosley 1993). 

It goes without saying that only the computation of a series of evolutionary tracks may 
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really provide a quantitative estimate of these additional important uncertainties. 

Before closing this section let us compare our final yields with the other ones available 
in the literature. We already published yields for 26 A1 in 2003 and 2004 (Limongi & Chieffi 
2003; Chieffi & Limongi 2004); those yields (the 26 A1 ones only), unfortunately, must be 
totally disregarded since, when we updated the full set of cross sections adopted in the 
evolutionary code, we found that the 26 Al(p, 7) 27 Si cross section present in the REACLIB 
database, and based on the Caughlan & Fowler (1988) formula, was wrong. Figure 5 shows 
the 26 A1 yields by WW95 (green filled squares), Rauscher et al. (2002) (cyan filled triangles), 
Thielemann et al. (1996) (magenta filled rhombs), Meynet et al. (1997) (red open dots), 
Palacios et al. (2005) (magenta open squares) and Langer et al. (1995) (green open rhombs) 
together with the present ones (marked as "this paper" - blue filled dots). The first thing 
worth noting is that our set of models is the only one that fully extends over the full range 
of the massive stars. How we have already said above, for historical reasons two parallel 
fields of research may be identified: the first one (mainly involving Woosley and coauthors), 
addressed the full evolution of stellar models (including the computation of the passage of 
the shock wave) but did not extend the explored mass interval above 35-4OM and a second 
one (that includes mainly European groups) involved in the study of the central H and He 
burning phases in a wide mass interval (up to 120M Q ) by including the mass loss and hence 
fully exploring the properties of the WR stars. The WW95 26 A1 yields show remarkable 
similarities with ours: both sets of yields show two maxima, a first one around 15-2OM and 
a second one around 3O-35M , the WW95 ones being slightly shifted towards more massive 
stars. We already discussed above how these two peaks are related to the temporal evolution 
of the C convective shell, and hence these strong similarities clearly point towards a quite 
similar presupernova evolution. The quite larger WW95 yields between 25 and 4OM are 
probably at least partially due to the contribution of the v process that we neglect. WW95 
quote the 26 A1 yield for the 25M without the v contribution and the value closely matches 
our yield. The Rauscher et al. (2002) yields are an upgrade of the WW95 for the solar 
metallicity; note that they do not show any more a strong peak around a 2OM , which 
means that something in the properties of the C convective shell has changed with respect 
to the older WW95 stellar models. Anyway, we feel confident to say that these three sets of 
yields lead to quite similar predictions. An additional set of yields for this mass interval was 
provided by Thielemann et al. (1996): these yields are much lower than all the other ones 
probably because they come from the evolution of just pure He cores computed by adopting 
a quite poor network in the hydrostatic burnings. Historically, stars more massive than, say, 
35M have been studied almost exclusively up to the end of the central He burning, so that 
only the 26 A1 present in the wind has been (widely) discussed in the literature. The latest 
26 Al yields in this mass interval are those provided by Langer et al. (1995), Meynet et al. 
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(1997) and Palacios et al. (2005). Note that, for homogeneity reason, only the non rotating 
models are shown. Figure 5 shows that the differences among these theoretical yields are 
confined within a factor of two. 



4. 60 Fe 

60 Fe is an unstable nucleus (its terrestrial half life is Tin — 1.5 x 10 6 y) that lies slightly 
out of the stability valley, its closest stable neighbor being 58 Fe. It is mainly produced by 
neutron capture on the unstable nucleus 59 Fe and destroyed by the (n, 7) process. Since 
its closest parent, 59 Fe, is unstable, the 59 Fe(n, 7) process must compete with the 59 Fe(/3~) 
decay to produce an appreciable amount of 60 Fe. An order of magnitude estimate of the 
neutron densities needed to cross the 59 Fe bottleneck may be derived by equating the (11,7) 
and f3~ decay rates: a value of the order of 3 x 10 10 11/cm 3 is obtained at temperatures 
lower than Log(T)=8.7. Above this temperature the half life of 59 Fe (whose terrestrial value 
is ~ 44 d) reduces systematically, dropping to ~ 6 d at Log(T)=9.0 and to ~ 1 hour at 
Log(T)=9.6, while the 59 Fe(n, 7) cross section varies very mildly with the temperature. As 
a consequence, the neutron density needed to cross the 59 Fe bottleneck steeply increases 
with the temperature: it raises to ~ 3 x 10 11 n/cm 3 at Log(T)=9 and to ~ 6 x 10 12 n/cm 3 
at Log(T)=9.3. Also the half life of 60 Fe has a steep dependence on the temperature: it 
remains of the order of the terrestrial value up to Log(T)=8.6 and then drops to half a year 
at Log(T)=9, to ~ 1.5 d at Log(T)=9.2 and to ~ 14 m at Log(T)=9.4. In spite of this strong 
dependence of its half life on the temperature, and contrarily to what happens to 26 A1, 60 Fe 
is mainly destroyed by the (11,7) reaction because the 60 Fe(n,7) rate always overcomes the 
decay rate at the extremely high neutron densities required to produce it. A temperature of 
the order of 2 billion degrees represents an upper temperature limit for the synthesis of 60 Fe 
because above this temperature the (7,11) and the (7, p) photodisintegrations of both 60 Fe 
and 59 Fe become tremendously efficient. Such an occurrence limits a possible 60 Fe production 
to the He, C and Ne burning phases. 

In He burning, where the main neutron donor is the 22 Ne(a,n) 25 Mg process, a temper- 
ature of the order of 4 x 10 8 K would be required to reach the threshold neutron density of 
3 x 10 10 n/cm 3 . In central He burning the temperature remains well below 3 x 10 8 K so that 
the neutron density never exceeds 10 7 n/cm 3 and no appreciable production of 60 Fe occurs. 
In shell He burning, vice versa, the temperature may raise up to and above 4 x 10 8 K so 
that even a large amount of 60 Fe may be synthesized. To understand its synthesis in this 
phase it is necessary to briefly remind a few firmly established characteristics of the shell 
He burning. At variance with the central H burning, where the mass size of the convective 
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core shrinks together with the H content, central He burning is characterized by a convective 
core that advances progressively in mass (or remains stable at most) as the He is depleted. 
As a consequence, the He profile left by the central He burning shows a sharp discontinuity 
in correspondence to the location of the maximum size of the convective core and the He 
convective shell forms outside this discontinuity in a region where the He abundance is flat 
and equal to the one left by the H burning; since all the advanced burning phases are very 
rapid with respect to the shell He burning timescale, only a modest amount of He is burnt 
in the shell before the core collapse so that its final abundance barely drops below 0.9 by 
mass fraction. The temperature at the base of this convective shell never raises enough to 
make the 59 Fe(n, 7) process competitive with respect to the 59 Fe(/3~) decay so that no 60 Fe 
is synthesized by these stars. Such a feature holds up to the first mass that becomes a WR 
star of the WNE/WCO kind (4OM ). In fact, these stars experience such a strong mass loss 
that they firstly lose all their H rich envelope and then they continue eroding the He core 
up to the moment of the core collapse. Since the properties of the He burning depend on 
the He core mass size, the core "feels" a progressively smaller mass as the He core reduces 
and hence the mass size of the convective core reduces accordingly: a He profile is left in 
these stars at the central He exhaustion and its shape will depend on the balance between 
the speed of the central burning and the mass loss rate. The He convective shell forms, in 
these stars, within the region of variable He abundance and hence the well known problem 
arises if the Schwarzschild criterion or the Ledoux one must be used to determine if (and 
on which timescale) a convective region forms. We do not want to address such a question 
here, mainly because a full non local treatment of the mixing is still missing and only semi 
phenomenological (and parametric) solutions have been proposed over the years to treat 
these layers. We decided, vice versa, to explore the two limiting cases: the first one in which 
the Schwarzschild criterion is adopted and a second one in which the Ledoux criterion is 
used. Columns 8 and 9 in Table 3 show the amount of 60 Fe that is produced in the two 
limiting cases: it is evident that the presence of a fully convective region largely increases 
the amount of 60 Fe synthesized by the shell He burning. Typical neutron densities reached 
in the He shell in these very massive stars range between 6 x 10 10 n/cm 3 and 1 x 10 12 n/cm 3 . 
A rather modest, but worth being mentioned, additional amount of 60 Fe is also produced in 
the radiative tail of the He burning shell of the less massive stars. These stars experience a 
significant contraction of the core (including the radiative tail of the He shell) after the end 
of the central Oxygen burning phase that raises the temperature (and hence the neutron 
density) enough to allow the synthesis of some 60 Fe (column 7 in Table 3). As the initial 
mass increases, this phenomenon disappears because the region external to the CO core be- 
comes progressively more insensitive to the evolution of the deep interior of the star (mainly 
because of the presence of a very active and stable C convective shell). The contribution of 
the He convective shell to the synthesis of 60 Fe is shown in Figure 6 as green filled rhombs 
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when the Schwarzschild criterion is adopted and as green empty rhombs when the Ledoux 
criterion is adopted. 

C burning behaves similarly to the He burning. In central C burning the neutron density 
does not exceed a few times 10 7 n/cm 3 , so that also in this case no 60 Fe may be produced. 
The main neutron donor is, once again, the 22 Ne(o;,n) 25 Mg process and the reason why 
it can't provide a high neutron rate is, in this case, the relatively low concentration of a 
particles provided by the 12 C( 12 C,o;) 20 Ne reaction. In analogy with the He burning, the 
larger burning temperature (T > 10 9 K) at which the shell C burning occurs, allows a much 
larger production of a particles that translates into a much higher neutron density and 
hence in a conspicuous 60 Fe production. The typical neutron density we find ranges between 
6 x 10 11 n/cm 3 and 2 x 10 12 n/cm 3 in the mass range 20 to 120 M while it drops significantly 
at lower masses. The mild dependence of the neutron density at the base of the C-convective 
shell on the initial mass (for M>17M ) explains why the 60 Fe produced by the C convective 
shell (blue filled triangles in Figure 6 and col. 10 in Table 3) increases only modestly with 
the initial mass. It must be noted, at this point, that the presence of a convective shell 
plays a crucial role also for the synthesis of 60 Fe. In fact it has the double responsibility of 
bringing new fuel (a particles and 22 Ne) in the region where the active burning occurs and 
simultaneously of bringing the freshly made 60 Fe outward in mass, i.e. at lower temperatures 
where the neutron density becomes negligible and its half life increases substantially. The 
30M Q model is the best example of the importance of the convective burning: in fact, the 
minimum in the yield of the 60 Fe produced by the C shell (see Figure 6) corresponding to the 
3OM model is the direct consequence of the lack of an efficient C convective shell lasting 
all along the advanced evolutionary phases (see the previous section). 

Ne burning would produce 60 Fe too, because of the large abundance of the a particles 
and neutrons provided in this case by the 23 Na(a,n) 26 Al reaction, but the lack of an extended 
and stable convective shell lasting up to the moment of the explosion prevents the build up 
of a significant amount of 60 Fe. In fact, the 60 Fe produced by the advancing Ne burning shell 
drops quickly to zero together with the local Ne abundance because either the production 
channel dries up and the f3~ decay speeds up consistently. By the way, the Ne convective 
shell that forms in the 14 and 15M after the end of the central Si burning, and that is 
responsible of the strong peak in the 26 A1 yield produced by these masses (see Figure 2), 
does not lead to a similar effect for the 60 Fe because the neutron density remains too low to 
cross the 59 Fe bottleneck. 

The last episode of synthesis of 60 Fe occurs when the blast wave crosses the mantle of 
the star on its way out, in the region where the peak temperature is of the order of 2.2 x 10 9 
K and hence roughly in the same region where the (explosive) synthesis of 26 A1 occurs. In 
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most cases such a temperature is reached either at the base or within the C convective shell 
so that the amount of explosive 60 Fe produced will depend on the local abundances of 20 Ne, 
12 C, 23 Na, 22 Ne left by the last C convective shell episode and on the final MR relation at 
the moment of the core collapse. The red filled squares in Figure 6 and col. 11 in Table 
3 show the contribution of the explosive burning to the synthesis of 60 Fe as a function of 
the initial mass. Two local minima corresponding to the two inversions in the MR relation 
are quite evident also in this case (see the previous section). In the mass interval 11 to 
15M Q , the explosive contribution to the synthesis of 60 Fe is almost negligible because of the 
low concentration of the nuclei that are necessary for its synthesis. The 13M constitutes, 
however, a striking exception because a large amount of 60 Fe is synthesized by the blast wave 
in this case. The reason is that the peak temperature of 2.2 x 10 9 K occurs beyond the outer 
border of the C convective shell where the abundances of 12 C and 22 Ne, in particular, are 
much higher than in the C convective shell. Since such a phenomenon occurs for just a single 
mass, one could question its reality. A close look at the models shows that, while the radius 
at which the blast wave synthesizes 60 Fe in these low mass stars remains roughly constant 
(~ 0.013R Q ), the radius of the outer border of the last C convective shell at the moment 
of the core collapse shows a non monotonic trend with the initial mass and, in particular, 
a minimum value of O.O12R for the 13M Q occurs. Such a minimum corresponds to a local 
minimum in the mass size of the last C convective shell as well. It is not easy to understand 
why the last C convective shell has a minimum mass size for this mass because the temporal 
evolution and mass (and radial) extension of the various convective episodes that follow each 
other depend on a non linear interplay among various factors (the main ones being the mass 
sizes of both the He and the CO core and the amount of C left by the central He burning). 
To better assess the range of masses for which this peculiar explosive production occurs, the 
computation of a finer mass grid would be required: at present we can just say that it is 
confined somewhere between 12 and 14M Q . Let us eventually note that no 60 Fe is lost by 
the stellar wind because, in the present computations, no fraction of the He convective shell 
is lost through the wind during the advanced WR phase. 

The black filled dots in Figure 6 and col. 6 in Table 3 show the total yield of 60 Fe as a 
function of the initial mass in the case in which the Schwarzschild criterion is assumed in the 
He convective shell. Its trend obviously reflects the trends of the main contributors at the 
various masses. There is a global direct scaling with the initial mass and a quite monotonic 
behavior: the two exceptions are the 13M and the 30M Q . Below the 4OM the total yield 
is dominated by the contribution of the C convective shell while above this mass it is the He 
convective shell to play the major role. The explosive burning almost always plays a minor 
role. The black open dots refer to the case in which the Ledoux criterion is assumed in the 
He convective shell: its contribution to the global synthesis of 60 Fe is strongly reduced in 
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this case. 

Similarly to the 26 A1 case, the treatment of the convective layers and the cross section 
of the 12 C(o;,7) 16 are strong sources of uncertainty in the prediction of the 60 Fe yield (see 
the previous section). The situation is even worse in this case because the cross sections of 
all the key processes involved in the synthesis of 60 Fe are purely theoretical: no experimental 
data exist for the 59 Fe(n,7) and the 60 Fe(n,7) rates as well as the dependence of their f3~ 
decay on the temperature. Let us eventually note that, since 60 Fe is mainly synthesized by 
the C and He shells, it may well be considered a pure secondary element because both the 
target (Fe) and the bullet (neutrons) depend on the initial metallicity (the neutrons, in fact, 
mainly come from the 22 Ne(a,n) process). Therefore also the rate of this process as well as 
the adopted initial [CNO/Fe] ratio (scaled solar or not) will significantly influence the final 
60 Fe outcome . 

60 Fe yields have been published, as far as we know, by Timmes et al. (1995), (they are 
part of the WW95 large database) and by Rauscher et al. (2002) and are shown in Figure 7 
as red filled squares and black filled triangles, respectively. Our yields obtained for the two 
extreme cases are shown as filled and open blue dots. For sake of completeness also the yields 
obtained by adopting the Langer (1989) mass loss rate are shown (see next section). The 
first thing worth noting is that our yields are the only ones available above 4OM . Second, 
a good overall agreement exists between the WW95 and our yields: both sets show a peak 
towards the lower end (13M in both cases) and a minimum (2OM in WW95 and 3OM in 
our models) that clearly indicate how the two sets of models behave similarly (at least 
qualitatively). As for Rauscher et al. (2002), their yields are systematically and significantly 
larger than both the WW95 and our yields. We do not have an definite explanation for such 
a difference since our models, computed with the latest input physics available today, show a 
remarkable agreement with the WW95 ones that have been computed with the input physics 
available at that time and not with the Rauscher et al. (2002) one that has probably been 
computed with a much closer input physics. Let us eventually remark that our 60 Fe yields 
already published in 2003 and 2004 (Limongi & Chieffi 2003; Chieffi & Limongi 2004) are in 
very close agreement with the present ones. 



5. Dependence of the yields on the Mass Loss rate adopted beyond the WNL 

phase 

In the previous sections we have discussed at some extent the yields of both the 26 A1 and 
60 Fe as a function of the initial mass as well as their dependence on some nuclear processes 
and some physical phenomena (like the efficiency of the convective mixing). In this section 
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we show how the yields of these two nuclei depend on the mass loss rate beyond the WNL 
phase. As it is well known, mass loss is still one of the main uncertain physical phenomenon 
that influence the evolution of a star. To explore the influence of the mass loss in the 
WNE+WCO phases we have recomputed the evolution of the stellar models that enter the 
WNE (i.e. stars in the range 40 to 120M Q ) by adopting the mass loss rate proposed by 
Langer (1989) and widely adopted in the past. This older mass loss rate is much stronger 
than the Nugis & Lamers (2000) one and it scales as M 2 ' 5 , so that all the very massive 
models tend now to converge towards a similar quite small final mass (Column 2 in Table 
5). Moreover, since such a strong reduction of the total mass occurs during the central He 
burning, all these models tend to have a quite similar CO core mass (column 3 in Table 
5) and hence quite similar advanced burning phases. This is the reason why the 26 A1 and 
60p e p ro d ucec i by the advanced burnings (hydrostatic plus explosive) tend to become quite 
insensitive to the initial mass (Columns 5 and 7 in Table 5). Of course the 26 A1 lost by the 
wind (Column 4 in Table 5) preserves, vice versa, a strong dependence on the initial mass 
because it is ejected in the earliest phases of the evolution of these stars. Note that the 
26 Al ejected by the wind is slightly larger in the Langer case because a larger fraction of 
the He core is ejected before being reached by the He burning. The total amount of 26 A1 
obtained by adopting the Langer (1989) mass loss rate (Column 4 in Table 5) is compared 
to the one obtained by adopting the Nugis & Lamers (2000) mass loss rate in Figure 8 (left 
axis). A similar comparison for the 60 Fe is shown in the same Figure (right axis). It is clear 
from these comparisons that while the 26 A1 yield is reduced by roughly a factor of two, the 
60 Fe one shows a dramatic dependence on the adopted mass loss rate. The reason is that 
the 60 Fe yield is dominated, in these very massive stellar models, by the He convective shell 
contribution (see the previous section) and the very strong mass loss implied by the Langer 
prescription completely kills such a contribution because it reduces enormously the final CO 
core mass. Hence, also the stability criterion adopted in the He convective shell does not 
play any significant role in this case. 



6. Comparison with the observations 

The clumpy and patchy distribution of the 26 A1 in the galactic plane shown by COMP- 
TEL, coupled to the strong correlation between its all-sky map and the 53 GHz free- free 
emission map, point towards a massive star parent for the 26 A1 and, more specifically, to- 
wards the subset of massive stars that are also responsible for the Lyman continuum photons 
that power the ionization of the interstellar medium. By comparing the 53 GHz free-free 
emission to the 26 A1 decay map, Knodlseder (1999) determined also that the scaling factor 
between the 1.8 MeV gamma ray line (7i.s) flux and the Lyman continuum photons (LCP) 
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flux remains roughly constant all over the galactic plane. He quantified this occurrence by 
identifying the average amount of 26 A1 (Y^ 7V = 1O _4 M ) that must be associated to the 
LCP flux (Q ) of an equivalent 07 V star (Log(Q ) = 49.05) in order to obtain the observed 
correlation between the two aforementioned all-sky maps; see Knodlseder (1999) for details. 
For reasons of convenience we prefer to reinterpret these data in terms of number of 71.8 
photons per ionizing photon (GPL); this number, that we name for convenience RgxL, re- 
sults to be equal to 1.25 x 10~ n GPL. Let us remind once again that the determination of 
PgxL does not require any explicit assumption about the sources of the ionizing photons and 
of the 26 A1 but only that the 26 A1 nuclei and the LCP share the same spatial distribution 
(Knbdlseder 1999). 

We have discussed in the previous sections essentially three sets of models that differ 
because of the mass loss rate adopted in the WNE+WCO phases and because of the stability 
criterion adopted in the He convective shell. For sake of clarity, the set of models computed 
by adopting the Nugis & Lamers (2000) mass loss rate in the WNE+WCO phases and the 
Schwarzschild criterion in the He convective shell will be referred to as the NLOO models, the 
one computed with the same mass loss rate but the Ledoux criterion in the He convective 
shell will be referred to as the NLOOL models while the set computed by adopting the Langer 
(1989) mass loss rate will be referred to as the LA89 models. Since all these different choices 
concern only stars more massive than 35M , all these sets share the same stellar models 
between 11 and 35M Q . 

From the theoretical side we can easily determine the RgxL factor (s) for our grid of 
models. Let us remind that at the base of this derivation there is the usual assumption of 
a steady state condition, i.e. of the equivalence between destruction and production rates: 
this implies that the observed decay rate of the 26 A1 is equal to the average ejection rate 
of this nucleus from the presently evolving stars. We have estimated the number of Lyman 
continuum photons produced by each of our models by means of table 3 of Schaerer & de 
Koter (1997); the two following relations provide both the average Log(Qo) and the lifetime 
of a star as a function of the initial mass (for our grid of models): 



Log 10 (Q ) = 34.4906 + 14.90772 x Log w (M) - 3.64592 x Log w (M) 2 (LCP) s _1 
Log 10 (t) = 9.598 - 2.879 x Log w (M) + 0.6679 x Log w {M) 2 (yr) 

Note that these relations are obviously identical for the three sets of models. Figure 9 
shows the trend of both Logio(Qo) (red filled dots, left axis) and the average 26 A1 production 
rate (first right axis) with the initial mass. The blue filled and open squares refer to the 
NLOO and the LA89 models, respectively. By the way, the average 26 A1 production rate is 
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simply the total amount of 26 A1 produced by a star divided by its lifetime. The RqxL factor 
per each stellar mass is also shown in the same Figure and its scale is on the second right 
axis. The black filled and open rhombs refer, respectively, to the NLOO and LA89 cases. RgxL 
remains roughly constant above 3OM because both Logio(Qo) an d the 26 A1 production rate 
run almost parallel; below this mass, on the contrary, the Lyman continuum photons drop 
down faster than the 26 A1 production rate, so that Rq x l increases significantly. Note that 
R-GxL varies by roughly a factor of two between the NLOO and the LA89 stellar models. 

The data shown in Figure 9 cannot be directly used to interpret the observed galactic 
value of RgxL because stars of different masses contribute simultaneously to both the LCP 
and the 71. § fluxes, and hence a proper initial mass function (IMF) must be taken into 
account. The IMF is usually expressed as a power law, namely: 

n{m) = K x m~( 1+x ^ 
where n(m) is the number of stars of mass m. The constant K is fixed by the condition: 

rM TO p 

N TO T = K x m- {l+x) dm = 1 

J Mbot 

where Mbot an d Mtop represent the limiting masses that bracket the range of the stars 
that explode as core collapse supernovae. As for the lowest mass that explodes as a core 
collapse supernova, we have adopted the smallest star in our grid, i.e. the 11M . One 
could question that even less massive stars do explode as core collapse supernovae, but this 
is practically uninfluent with respect to the computation of the RgxL parameter because 
below 11M both the LCP and the 26 A1 drop practically to zero. On the other hand, the 
choice of the upper limit may play a role. Figure 10 shows the theoretical RgxL for a range 
of possible values of the power law slope, x, and three different values of the upper mass 
limit Mtop- The blue solid, red dotted and green dashed lines refer to three different choices 
of the upper mass limit: 120, 60 and 4OM respectively. The thick and thin lines refer to 
the NLOO and LA89 models respectively. It is quite evident that, fortunately, RgxL does not 
depend significantly on the adopted mass loss rate and also the dependence on Mtop is quite 
modest; in fact, even a drastic change of Mtop from, e.g., 120 to 4OM affects RgxL by no 
more than a factor of two. The reason for such a low dependence of RgxL on Mtop is that 
above 3OM the ratio between the 26 A1 production rate and Lyman continuum luminosity 
remains roughly constant (see Figure 9) for all sets of models. 

The horizontal thick solid black line marks the observed galactic value derived by 
Knodlseder (1999), but expressed in terms of RgxL, while the shaded area reflects the range 
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of values that correspond to an uncertainty of a factor of two in the observed value. Kroupa 
& Weidner (2003) and Kroupa (2004) have recently rediscussed the IMF slope for massive 
stars and came to the conclusion that it must be larger than, at least, x=1.8. If we take this 
specific value (marked by the thin vertical black dashed line in Figure 10) as a conservative 
representation of the actual IMF slope of these stars, it turns out that a good fit to the 
galactic RqxL exists for any Mtop mass ranging between, at least, 120 and 4OM and both 
sets of models. The fit improves for steeper IMF slopes. We therefore conclude that either 
the NL00 and the LA89 sets of models can naturally explain and reproduce the constant 
galactic RgxL- 

The good fit to this parameter, however, does not carry any information, by itself, about 
the total amount of 26 A1 present in our Galaxy. In order to determine such an amount by 
means of RgxL 5 it is necessary to know the total Lyman continuum luminosity provided by 
our Galaxy (Qmw)- Once this value is known, the global 26 A1 production rate, P M w> is given 
by ?mw = Qmw x R-GxL and the present steady state abundance of 26 A1 in our Galaxy would 

i • i t-i 26 A1 26 A1 

be given by P M $ x r dcc A1 . 



Recent determinations of the galactic Lyman continuum luminosity range between 
Qmw = 3.5 x 10 53 photons s _1 (Bennett et al. 1994) and Qmw = 2.6 x 10 53 photons s _1 
(McKee & Williams 1997) with about 50% uncertainty in both cases (see the cited papers for 
details). If one would adopt the highest value (Bennett et al. 1994), an IMF slope x=1.8 and 
Mtop — 120M Q , the present steady state abundance of 26 A1 in our Galaxy would amount to 
1.97M Q for the NL00 models and to 1.70M Q for the LA89 case. Figures 10 (right scale) and 11 
show the trends of the 26 A1 abundance, of the SFR, of the number of Type II supernovae and 
of the Type Ibc supernovae as a function of the slope x for a Qmw = 3.5 x 10 53 photons s _1 . 
In Figure 11, the thick lines refer to an Mtop — 120M Q while the thin lines refer to an 
Mtop — 6OM . Since all these relations scale linearly with Qmwj the adoption of the value 
provided by (McKee & Williams 1997) would simply imply a downward scaling by a factor 
equal to 2.6/3.5. The quantities shown in Figure 11 obviously do not depend at all neither on 
the adopted mass loss beyond the WNL phase nor on the stability criterion adopted in the 
He convective shell. Figure 12 shows the total amount of 26 A1 ejected by a stellar generation 
(solid lines - right axis) together with the percentage contributions of the Type II supernovae 
(dot dashed lines), the Type Ibc supernovae (dashed lines) and the WR stars (dotted lines), 
all as a function of the IMF slope x: the thick and thin lines refer to the NL00 and the LA89 
cases, respectively. The observed trends may be readily understood by reminding that a) 
the LA89 models predict sistematically less 26 A1 than the NL00 ones (see Table 5) because 
of the reduced production both in the C convective shell and Ne explosive burnings (see 
the previous section) and b) the total amount of 26 A1 ejected by the Type II supernovae is 
identical in both sets because they share the same stellar models below 4OM . Since the 
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integration over the IMF may depend on the adopted upper mass limit (Mtop), we show 
in Figures 13 (for the NLOO case) and 14 (for the LA89 case) the same quantities shown in 
Figure 12 as a function of Mtop and two specific IMF slopes, i.e. x=1.8 and x=2.5. Both 
these Figures clearly show that, indipendently on Mtop an d on the mass loss rate adopted 
beyond the WNL stage, the largest contribution to the total 26 A1 produced by a stellar gen- 
eration comes from the Type II supernovae, while both the WR and the Type Ibc contribute 
at most (complessively) for a 30% in the most favoreable case. The present result, i.e. that 
most of the 26 A1 comes from the Type II supernovae must be interpreted as a property of 
our sets of models. We have already shown in Figure 5 that the yields produced by different 
authors differ even significantly in the mass range they share; it is therefore very probable 
that they could have been different also in the mass interval not explored by other groups so 
that even the relative contributions of the various mass intervals could change significantly 
from one author to another. In other words, it would be a mistake to use a set of models 
that extend only up to 30- 40M Q to infer the 26 A1 produced by a generation of massive stars 
basing such a choice on the properties of the present models. It is not known a priori the role 
that would play the more massive stars within other sets of models. In this respect it is of 
overwhelming importance that also other groups compute grid of models extended in mass 
at least as much as in the present ones, because only a comparison between indipendent 
computations would allow a better understanding of the role really played by the various 
mass intervals to the global budget of the 26 A1. It goes without saying that the same holds 
for the 60 Fe. 

Both RHESSI and SPI (INTEGRAL) have reported a measurement of the 60 Fe/ 26 Al 
flux ratio towards the central radiant of our Galaxy (see the introduction). The latest 
values, at the time of writing, are (0.17 ± 0.05) (RHESSI, Smith 2005) and (0.11 ± 0.03) 
(SPI, Harris et al. 2005). Since the two experimental data show a compatibility range around 
60 Fe/ 26 Al=0.14, we decided to tentatively (and arbitrarily) adopt this value as representative 
of both the experiments and to consider as a typical error the semi difference between the 
quoted values. This representative value is shown as a horizontal black thick line in Figure 
15, while the two experimental values are shown as thin black solid lines in the same Figure. 
Figure 15 also shows our theoretical predictions for the three sets of models, namely the 
NLOO (red solid line), the NL00L (blue dashed line) and the LA89 (green dotted line) and 
M T op=12OM . This Figure shows that the NLOO models predict a too large 60 Fe/ 26 Al 
flux ratio for any slope in the explored range. Such a large flux ratio is due to the very 
large amount of 60 Fe produced in the He convective shell in stars more massive than 35M Q . 
The adoption of the Ledoux criterion in the He convective shell significantly dumps out the 
60p e production in this region so that the theoretical flux ratio reduces somewhat though it 
remains larger than the observed value. The fit improves as the IMF slope increases because 
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the contribution of the more massive stars reduces correspondingly. The LA89 models, vice 
versa, predict an 60 Fe/ 26 Al flux ratio in excellent agreement with the observed value over 
a very large range of IMF slopes, and this is the consequence of the very strong mass loss 
that prevents either the production of 60 Fe in the He convective shell and reduces also its 
production in the C convective shell and in the Ne explosive burning. Figure 16 shows, for 
two selected IMF slopes, the dependence of this flux ratio on Mtop- Only the NLOO models 
show a strong dependence of the flux ratio on M TO p and this is due to the fact that this 
is the only set for which a substantial contribution to the 60 Fe comes from the stars more 
massive than 4OM . 

A paper devoted to the presentation of a new generation of 26 A1 yields cannot escape 
the temptation to reanalyze the fit to 7 2 Velorum, the closest WR star (WR11) presently 
known. As far as we know, the latest very detailed and comprehensive analysis of this WR 
star has been presented by Oberlack et al. (2000) and we refer the reader to this paper 
and to the references therein for a detailed discussion of both the observational aspects 
and the theoretical problems involved in the interpretation of this star (either in the single 
and binary stars scenarios). Here it suffices to say that the typical initial mass quoted for 
this star is 60 ± 15M , that its initial metallicity was very probably solar and the binary 
system is wide enough that probably the evolution of 7 2 Velorum was not affected much by 
the binarity environment. No 26 A1 has been detected around this star but an upper limit 
of 6.3l L4 x 1O _5 M has been derived by Oberlack et al. (2000). This value represents a 
problem because the theoretical models predict a much larger amount of 26 A1. Meynet et al. 
(1997) find that a 60M Q solar metallicity star ejects 1.49 x 1O~ 4 M of 26 A1 through the wind 
while Palacios et al. (2005) quote a value of 1.3 x 1O~ 4 M for the same initial mass; also the 
older models by Langer et al. (1995) predict an amount of 26 A1 equal to 1.2 x 1O~ 4 M for 
a 5OM star. Also for a mass as low as 4OM the models predict rather large values that 
are only barely compatible with the quoted upper limit: 5.5 x 1O~ 5 M (Meynet et al. 1997) 
and 5.09 x 1O~ 5 M (Langer et al. 1995). Such a failure in reproducing the closest WR star 
is very embarrassing because this is a member of the class of stars that are very probably 
responsible for the production of the bulk of the 26 A1 presently in the disk of our Galaxy. 
An inspection to Table 3 shows that the wind component of the 26 A1 yield in our models 
is roughly a factor of two lower than that obtained by the other authors, so that a large 
part of the discrepancy is removed. In addition to this, there is another important point 
that has been neglected up to now: the timescale over which the 26 A1 is injected into the 
interstellar medium. Figure 17 shows the cumulative abundance of 26 A1 present in the ejecta 
as a function of the age of the star and the four panels refer to the 40, 60, 80 and 12OM : 
the blue dashed thick lines show the case in which the 26 A1 is assumed to accumulate in the 
interstellar medium without decaying (unrealistic, but usually adopted, ideal case in which 
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the 26 A1 is assumed to start decaying after it has been completely ejected) while the red solid 
thick lines represent the real total amount of 26 A1 one would find in the interstellar medium 
as a function of time (taking into account the fact that the 26 A1 starts to decay as soon as 
it is ejected in the interstellar medium). The gray area mark the temporal phase in which 
each star appears as a WCO star. If 7 2 Velorum is now just at the beginning of its WCO 
phase (the worst case because the real 26 A1 abundance is at its maximum) the predicted 
amount of 26 A1 would be 15 to 20% smaller than the value shown in Table 3 for masses in 
the range 40 to 120M Q . In particular, an abundance of 5.82 x 1O _5 M is predicted in the 
interstellar medium in the case of a 60M Q model, a value well below the upper limit quoted 
at present. This means that at least the full range of masses between 4OM and 6OM is 
now compatible with the observed upper limit. There is even room for an initial mass larger 
than 6OM . We can't be more precise at the moment because our nearest grid mass is the 
8OM . The adoption of the LA89 models would not change these results appreciably. Hence 
we conclude that a proper treatment of the injection of the 26 A1 ejected into the interstellar 
medium coupled to our new yields removes the longstanding discrepancy between the upper 
limit quoted for 7 2 Velorum and the theoretical predictions. 



7. Conclusions 

We have extensively discussed the production sites of the two gamma ray emitters 26 A1 
and 60 Fe over the range of massive stars that may contribute significantly to either the 
production of these nuclei and to the Lyman continuum luminosity, i.e. the range 11 to 
12OM . These theoretical predictions fully cover all the evolutionary phases from the pre 
main sequence to the explosive burnings. At variance with current ideas, 26 A1 is mainly 
produced by the Ne/C explosive burning over most of the mass interval analyzed. The main 
production site of the 60 Fe, vice versa, strongly depends on the adopted mass loss rate. In 
the LA89 case the main 60 Fe producer is always the C convective shell while in the NL00 case 
the main producer is still the C convective shell for masses lower than 4OM while above 
this mass a strong contribution comes from the He convective shell. Since the He convective 
shell forms, in these stars, in a region where a gradient in the He abundance is present, its 
contribution to the synthesys of 60 Fe depends significantly on the adopted stability criterion 
(Schwarzschild or Ledoux). We have used these yields to address the problem of the diffuse 
abundances of 26 A1 and 60 Fe in the Galaxy. 

The discovery of an almost constant flux ratio all over the galactic disk between the 
Lyman continuum photons (derived from the 53 GHz free-free emission all-sky map) and the 
7i.8 MeV line (that we express as RqxL, number of 71.8 MeV photons Per Lyman continuum 
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photon -GPL), plus the recently determined 60 Fe/ 26 Al flux ratio towards the inner radiant 
of the Galaxy, constitute in our opinion the two key experimental data a set of models 
must satisfy. The reason is that, once they are assumed to be equally spatially distributed, 
their ratio is independent on their spatial distribution. In addition to that, also the highly 
uncertain SFR fortunately does not play any role in these ratios (once a steady state is 
assumed). The experimental value of the RgxL parameter has been determined by Knodlseder 
(1999) (even if in different units) to be Log(RcxL) = —10.9 ± 0.3 while the value of the 
60 Fe/ 26 Al flux ratio (averaged between the RHESSI and SPI values) is of the order of 0.14 ± 
0.03. As for the first ratio, RgxL, its constancy along the galactic longitude could hide 
additional information. Different longitudes probably map combinations of different stellar 
populations (the smallest longitudes probably mapping the region with the highest average 
metallicity): hence the constancy of RgxL could imply a low dependence of this ratio on the 
initial metallicity. The same argument does not hold for the 60 Fe/ 26 Al flux ratio because in 
this case we just have one integrated measure towards the central radiant of the Galaxy. 

We have shown that the observed RgxL may be reproduced fairly well for a variety of 
slopes of the IMF and both sets of models, namely the NL00 and the LA89 ones. If we focus 
on a specific slope, e.g. x=1.8 (a value widely accepted to be reasonable for these massive 
stars) a maximum discrepancy of 0.2 dex is obtained, independently on the adopted M TO p 
and mass loss rate beyond the WNL phase. The fit even improves for steeper slopes. Within 
our sets of models, 26 A1 is mainly produced by stars less massive than, say, 35M , the total 
contribution of the more massive stars (through the wind of the WR stars plus the Type Ibc 
supernovae) never exceeding 30% of the total. Also this result depends only moderately on 
the adopted mass loss rate beyond the WNL phase and on Mtop- 

The theoretical 60 Fe yield, on the contrary, strongly depends on the mass loss rate 
adopted in the WNE+WCO phases mainly because mass loss controls the amount of 60 Fe 
produced in the He convective shell of the stars more massive than 35M during the very 
latest phases that precede their final collapse. Depending on the adopted mass loss rate, 
the predicted 60 Fe/ 26 Al flux ratio may range between 0.12 and 0.25 for a slope x=1.8. The 
models that better reproduce the observed ratio (~ 0.14) are the LA89 ones since the very 
strong mass loss rate proposed by Langer (1989) completely dumps out the contribution of 
the He convective shell and also reduces the contribution of the C convective shell because 
of the much smaller CO core mass. The dependence of the 60 Fe/ 26 Al flux ratio on both the 
IMF slope and Mtop is negligible in this case because of the overall modest contribution of 
the stars more massive than 35M to the global budget of the 60 Fe. The worst fit is obtained 
for the NL00 models since in this case either a strong contribution to the total yield comes 
from the He convective shell and also the contribution from the C convective shell is larger 
because the mass size of the CO core is definitely much large in this case. A steep IMF plus 
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an Mtop not larger than, say, 60M Q could anyway reconcile also this set of models with the 
observed flux ratio. 

Of course many other unceretainties affect both yields and we have shown in this paper 
how the 26 A1 yield depends, e.g., on several cross sections, initial abundances and size of 
the H-convective core. Also a distortion of the IMF towards the highest masses could lead, 
in practice, to a steepening of the slope. It must also be noted that all the present results 
have been obtained with solar metallicity stars while very probably at least the inner Galaxy 
has an average metallicity larger than solar. We actually do not know how the metallicity 
influences the yields of the 26 A1 and the 60 Fe and hence we leave such an exploration to 
another paper. The dependence of the 26 A1 wind component on the initial metallicity has 
been already discussed in a number of papers (Prantzos 1991; Meynet & Arnould 1993; 
Meynet et al. 1997) but unfortunately such a dependence is not representative of the global 
trend because we have shown in this paper that the 26 A1 present in the wind always represents 
a modest contribution to the total budget of the 26 A1. 

The determination of the total amount of 26 A1 present in the Galaxy requires the knowl- 
edge of an additional quite uncertain parameter, i.e. the SFR. Such a quantity is usually 
determined by means of either the number of core collapse supernovae per unit time or the 
total Lyman continuum luminosity of the Galaxy. Since it would be desirable to have a SFR 
linked as strongly as possible to the stars that mainly contribute to the Lyman continuum 
luminosity as well as to the 26 A1 and 60 Fe synthesys, the adoption of the total Lyman con- 
tinuum luminosity should be preferred. The core collapse supernovae, in fact, are dominated 
(in number) by the Type II supernovae whose range extends down to masses that do not 
contribute significantly to neither the Lyman continuum luminosity nor the 26 A1 nor the 
60 Fe. For example, the mass range 11 to 13M produces roughly 30% of the total number of 
Type II supernovae while their contribution to the synthesis of the 26 A1 is of the order of 6%; 
the situation would even worsen if the lower mass limit of the massive stars would reduce to 
10 or even 9M because these stars would strongly affect the frequency of the Type II su- 
pernovae without contributing significantly neither to the Lyman continuum luminosity nor 
to the 26 A1 or to the 60 Fe budgets. By normalizing to the total Lyman continuum luminosity 
estimated by Bennett et al. (1994) and an IMF slope x=1.8, we predict roughly 2M Q of 26 A1 
presently in the Galaxy in the NL00 case. The adoption of the LA89 models would reduce 
such an amount to 1.7M . 

We have also fitted 7 2 Velorum, the closest WR star. In particular we predict an 
abundance of 5.82 x 1O _5 M of 26 A1 (for a 6OM ), value compatible with the current upper 
limit of 6.31:1a x 1O~ 5 M . 

We did not address the fit to specific OB associations in this paper because it would 
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require an evolutionary synthetic code, as it has been correctly addressed by, e.g., Cervino 
et al. (2000), that is not available at present. 
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Fig. 1. — Temporal evolution of the central 26 A1 mass fraction during core H burning. The 
black solid, blue long dashed, red dotted and green short dashed lines refer to the 15, 30, 60 
and 120M Q models, respectively. The cyan dot-dashed line refers to a pure 26 A1 decay. 
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axis, green dotted lines) supernovae per century as a function of the IMF slope x for two 
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120M Q (lower right). The blue dashed thick lines show the case in which the 26 A1 is assumed 
to accumulate in the interstellar medium without decaying; the red solid thick lines represent 
the real total amount of 26 A1 one would find in the interstellar medium as a function of time, 
taking into account that the 26 A1 starts to decay as soon as it is ejected into the interstellar 
medium. The gray area mark the temporal phase in which each star appears as a WCO star. 
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H3 


CA88 




CA88 


NACR 
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RATH 


He3 


FKTH 






NACR 


NACR 




FKTH 


CA88 














He4 










NACR 


FKTH 


















Li6 
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Li7 
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BK 
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Be9 


NACR 
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NACR 


FKTH 
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BelO 
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FKTH 
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FKTH 














RATH 


BIO 
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FKTH 


CA88 


FKTH 




NACR 












Bll 


NACR 








FKTH 


FKTH 


FKTH 
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C12 


NACR 






KUNZ 


NACR 


CA88 
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C13 


NACR 


FKTH 
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FKTH 
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C14 
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CA88 
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RATH 
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FKTH 


RATH 




NACR 
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RATH 
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NACR 
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NACR 


NACR 


RATH 


BK 


FKTH 


NACR 


NACR 


KUNZ 
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NACR 
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FKTH 


CA88 


FKTH 


NACR 


FKTH 




FKTH 
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OD94 
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NACR 
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FKTH 


NACR 




NACR 


BK 




FKTH 




CA88 


FKTH 




OD94 
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FKTH 


FKTH 






FKTH 
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OD94 
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FKTH 


CA88 






NACR 


RATH 


FKTH 


FKTH 
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OD94 
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FKTH 


FKTH 
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OD94 
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FKTH 
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Ne20 
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RATH 
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NACR 


NACR 
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OD94 
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IL01 




RATH 


CA88 


RATH 


NACR 


BKRT 




NACR 




CA88 


BKRT 
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Ne22 


IL01 


CA88 
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RATH 


JAEG 
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Ne23 


RATH 
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BKRT 




OD94 
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RATH 


NACR 




RT 


RATH 
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OD94 




Na22 
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OD94 
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RT 


RATH 
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RT 


FFN8 




Sc43 
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RT 


RATH 
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RATH 


RATH 


RT 
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Sc44 
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FFN8 


FFN8 
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LPOO 
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RT 
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LPOO 


Sc48 
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RT 
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Sc49 
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BKRT 


LPOO 
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BKRT 
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RATH 


RATH 


RATH 


BKRT 


LPOO 


LPOO 


Ti50 


RATH 


RATH 


RATH 


RATH 




RATH 


BKRT 




RATH 
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RATH 


BKRT 




LPOO 


Ti51 


RATH 


RATH 


RATH 


RATH 




RATH 






RATH 




RATH 


BKRT 




LPOO 


V45 










RATH 




RT 


RATH 


RATH 


FKTH 


RATH 




FFN8 




V46 








RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 




V47 


RATH 
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RATH 


RATH 


RATH 


RT 


RATH 


RATH 
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RATH 
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LPOO 
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V52 
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RATH 


RATH 
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BKRT 




LPOO 


Cr48 




RATH 




RATH 
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RT 


RATH 


RATH 
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LPOO 
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LPOO 




Cr50 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Cr51 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 
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RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 


RATH 


RATH 


BKRT 
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LPOO 


Cr55 
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RATH 


BKRT 




LPOO 


Mn50 








RATH 


RATH 




RT 


RATH 


RATH 


RATH 
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LPOO 




Mn51 


RATH 
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RATH 


RATH 


RATH 


RT 
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LPOO 
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RT 
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BKRT 
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RATH 
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RATH 


RATH 


RATH 


RATH 




RATH 




RATH 


RATH 


RT 


LPOO 


LPOO 


Ni56 


FKTH 






RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 




LPOO 




Ni57 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


FFN8 


Ni58 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Ni59 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 


LPOO 


LPOO 


Ni60 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 


LPOO 


LPOO 


Ni61 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 


LPOO 


LPOO 


Ni62 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 


RATH 


RATH 


BKRT 




LPOO 


Ni63 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




RATH 


BKRT 




LPOO 


Ni64 


RATH 


RATH 


RATH 


RATH 




RATH 


BKRT 




RATH 




RATH 


BKRT 




LPOO 


Ni65 


RATH 




RATH 


RATH 




RATH 










RATH 


BKRT 




LPOO 


Cu57 










RATH 




RT 


RATH 


RATH 


FKTH 






FFN8 




Cu58 








RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 




Cu59 


RATH 


RATH 




RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 




Cu60 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Cu61 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Cu62 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Cu63 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Cu64 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


BKRT 


LPOO 


LPOO 


Cu65 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 




RATH 


RATH 


RT 


LPOO 


LPOO 


Cu66 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 








RATH 




BKRT 




RATH 


Zn60 




RATH 




RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 




LPOO 





en 















Table 1 


. — Continued 














isotope 


(P.7) 


(P>") 


(P.n) 


(a,7) 


(<*>p) 


(a,n) 


0,7) 


(n, P ) 


(n,a) 


(7>P) 


(7>«) 


(7,n) 


(/3 + ) 


W~) 


Zn61 


RATH 


RATH 




RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 




Zn62 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Zn63 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Zn64 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Zn65 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 


LPOO 


LPOO 


Zn66 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 






Zn67 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 


RATH 


RATH 


BKRT 






Zn68 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




RATH 


BKRT 






Zn69 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 








RATH 


BKRT 




RATH 


Zn70 


RATH 




RATH 


RATH 




RATH 


BKRT 










RT 






Zn71 


RATH 




RATH 


RATH 




RATH 












BKRT 




RATH 


Ga62 










RATH 




RT 


RATH 


RATH 


RATH 


RATH 




LPOO 




Ga63 


RATH 


RATH 






RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 




Ga64 


RATH 


RATH 


RATH 




RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Ga65 


RATH 


RATH 


RATH 




RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 


LPOO 


Ga66 


RATH 


RATH 


RATH 




RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




Ga67 


RATH 


RATH 


RATH 


RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




Ga68 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




Ga69 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 






Ga70 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




RATH 


RATH 


BKRT 




RATH 


Ga71 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 




RATH 




RT 






Ga72 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 








RATH 




BKRT 




RATH 


Ge64 














RT 


RATH 


RATH 


RATH 


RATH 




LPOO 




Ge65 




RATH 










RT 


RATH 


RATH 


RATH 


RATH 


RT 


LPOO 




Ge66 




RATH 










RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




Ge67 




RATH 










RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




Ge68 




RATH 






RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




Ge69 




RATH 






RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




Ge70 


RATH 


RATH 




RATH 


RATH 




BKRT 


RATH 


RATH 


RATH 


RATH 


RT 






Ge71 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 




Ge72 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 






Ge73 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 


RATH 


RATH 


BKRT 






Ge74 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




RATH 


BKRT 






Ge75 


RATH 


RATH 


RATH 


RATH 




RATH 


RT 








RATH 


BKRT 




RATH 


Ge76 


RATH 




RATH 


RATH 




RATH 


BKRT 










RT 






















Table 1 — Continued 
















isotope 


(P.7) 


(P: Q ) 


(P.n) 


0,7) 


( Q .P) 


(a,n) 


0,7) 


(n,p) 


(n,a) 


(7,P) 


(7,«) 


(7,n) 


09+) 


OS" ) 




Gc77 






RATH 


RATH 




RATH 












BKRT 




RATH 




As71 




RATH 




RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 




RATH 






As72 




RATH 




RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






As73 


RATH 


RATH 




RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






As74 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 




As75 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 








As76 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




RATH 


RATH 


BKRT 




RATH 




As77 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 




RATH 




RATH 




RT 




RATH 




Se74 


RATH 


RATH 




RATH 


RATH 




BKRT 


RATH 


RATH 


RATH 


RATH 










Se75 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 






Se76 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 








Se77 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 








Se78 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 


RATH 


RATH 


BKRT 








Se79 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




RATH 


BKRT 




TAKA 




Se80 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




RATH 


BKRT 








Se81 


RATH 




RATH 


RATH 




RATH 


RT 








RATH 


BKRT 




RATH 


00 


Se82 


RATH 




RATH 


RATH 




RATH 


BKRT 










RT 






Se83 






RATH 


RATH 




RATH 












BKRT 




RATH 




Br75 








RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 




RATH 






Br76 








RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Br77 


RATH 


RATH 




RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Br78 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Br79 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 








Br80 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 




Br81 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 




RATH 


RATH 


RT 








Br82 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 




RATH 




BKRT 




RATH 




Br83 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 




RATH 




RATH 




RT 




RATH 




Kr78 


RATH 


RATH 






RATH 




BKRT 


RATH 


RATH 


RATH 


RATH 










Kr79 


RATH 


RATH 


RATH 




RATH 




BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 






Kr80 


RATH 


RATH 


RATH 


RATH 


RATH 




BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 








Kr81 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 


TAKA 






Kr82 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 








Kr83 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 








Kr84 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 


RATH 


RATH 


BKRT 








Kr85 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




RATH 


BKRT 




TAKA 

















Table 1 — Continued 
















isotope 


(P.7) 


(P.«) 


(P.n) 


(«>7) 


(<*>p) 


(a,n) 


(n.7) 


(n, P ) 


(n,a) 


(7,P) 


(7,a) 


(7,n) 


(/3+) 


09") 




Kr86 


RATH 


RATH 


RATH 


RATH 




RATH 


BKRT 




RATH 




RATH 


BKRT 








Kr87 


RATH 




RATH 


RATH 




RATH 










RATH 


BKRT 




RATH 




Rb79 














RT 


RATH 


RATH 


RATH 


RATH 




RATH 






Rb80 














RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Rb81 




RATH 




RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Rb82 




RATH 




RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Rb83 


RATH 


RATH 




RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Rb84 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 




Rb85 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 








Rb86 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 




TAKA 




Rb87 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 




RATH 


RATH 


BKRT 




RATH 




Rb88 


RATH 


RATH 


RATH 




RATH 


RATH 








RATH 




BKRT 




RATH 




Sr84 


RATH 


RATH 






RATH 




BKRT 


RATH 


RATH 


RATH 


RATH 










Sr85 


RATH 


RATH 


RATH 




RATH 




RT 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 






Sr86 


RATH 


RATH 


RATH 


RATH 


RATH 




BKRT 


RATH 


RATH 


RATH 


RATH 


RT 








Sr87 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 








Sr88 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 






Sr89 


RATH 


RATH 


RATH 


RATH 




RATH 


BKRT 




RATH 


RATH 


RATH 


BKRT 




RATH 




Sr90 


RATH 


RATH 


RATH 


RATH 




RATH 


RT 




RATH 




RATH 


BKRT 




RATH 




Sr91 




RATH 


RATH 


RATH 




RATH 










RATH 


RT 




RATH 




Y85 














RT 


RATH 


RATH 


RATH 


RATH 




RATH 






Y86 














RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Y87 




RATH 




RATH 


RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Y88 




RATH 




RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 






Y89 


RATH 


RATH 




RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 








Y90 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




Y91 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 




RATH 


RATH 


RATH 


RATH 


RT 




RATH 




Zr90 


RATH 


RATH 




RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 










Zr91 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 








Zr92 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 


RATH 


RATH 


BKRT 








Zr93 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




RATH 


BKRT 




RATH 




Zr94 


RATH 


RATH 


RATH 


RATH 


RATH 


RATH 


BKRT 




RATH 




RATH 


BKRT 








Zr95 


RATH 




RATH 






RATH 


BKRT 








RATH 


BKRT 




TAKA 




Zr96 


RATH 




RATH 








BKRT 










BKRT 




RATH 




Zr97 






RATH 


















BKRT 




RATH 





Table 1 — Continued 



isotope (p,7) (p,a) (p,n) (0,7) (a,p) (a,n) (11,7) (n,p) (n,a) (7,p) (7,0) (7,11) (/?+) (/3 ) 



Nb91 


RATH 






Nb92 


RATH 




RATH 


Nb93 


RATH 


RATH 


RATH 


Nb94 


RATH 


RATH 


RATH 


Nb95 


RATH 


RATH 


RATH 


Nb96 


RATH 


RATH 


RATH 


Nb97 
Mo92 
Mo93 
Mo94 
Mo95 
Mo96 
Mo97 
Mo98 


RATH 


RATH 

RATH 
RATH 
RATH 
RATH 
RATH 


RATH 



RATH 
RATH 
RATH 
RATH 
RATH 



RT 


RATH 


RATH 


RATH 


RATH 




RATH 




RT 


RATH 


RATH 


RATH 


RATH 


RT 


RATH 


RATH 


BKRT 


RATH 


RATH 


RATH 


RATH 


RT 






BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 




TAKA 


RT 


RATH 




RATH 


RATH 


BKRT 




RATH 


RT 


RATH 




RATH 




RT 




RATH 




RATH 




RATH 




RT 




RATH 


BKRT 


RATH 




RATH 










RT 


RATH 


RATH 


RATH 




BKRT 


TAKA 




BKRT 


RATH 


RATH 


RATH 


RATH 


RT 






BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 






BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 






BKRT 


RATH 


RATH 


RATH 


RATH 


BKRT 










RATH 


RATH 


RATH 


BKRT 







Note. — BK = Bao et al. (2000). BKRT = Bao et al. (2000) in the energy range 5-100 keV; Rausher & Thielemann (2000) above this limit 
but rescaled to match the experimental values of Bao et al. (2000) at 100 keV. CA88 = Caughlan & Fowler (1988). FKTH = Thielemann et 
al. (1995), Reaction Rate Library REACLIB in which the experimental values are preferred whenever available, (for the experimental rates the 
references are reported in the library). IL01 = Iliadis et al. (2001). JAEG = Jaeger et al. (2001), recommended rate. KUNZ = Kunz ct al. 
(2002), adopted rate. LUNA = Formicola et al. (2004), LUNA (Laboratory, for Underground Nuclear Astrophysics) collaboration. NACR = 
Angulo ct al. (1999), NACRE. RATH = last version of the REACLIB provided by T. Rauscher and F. Thielemann and also available at the 
web site http://quasar.physik.unibas.ch/ tommy. RT = Rauscher & Thielemann (2000). OD94 = Oda ct al. (1994). FFN8 = Fuller, Fowler & 
Newman (1982,1985). LP00 = Langankc & Martinez Pincdo (2000). TAKA = Takahashi & Yokoi (1987) 



O 



Table 2. Basic evolutionary properties of the stellar models 



M ini 


M cnd 


t H 


tHc 


^t CX pl_ 


M Hc 


M CO 


Mbot&top(C sn ) 


vccn 


to 


t\VR 


t\VNL 


twNE 


twco 


Mq 


M 


Myr 


Myr 


yr 


M 


M 


M 


mass fraction 


yr 


yr 


yr 


yr 


yr 


11 


10.56 


20.11 


1.55 


7.32E-4 


3.473 


1.749 


1.5-1.739 


0.359 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


12 


11.49 


17.44 


1.30 


6.24E-4 


3.912 


1.971 


1.6-1.916 


0.361 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


13 


12.01 


15.40 


1.14 


4.58E-4 


4.363 


2.212 


1.8-2.011 


0.350 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


14 


12.77 


13.81 


1.01 


3.73E-4 


4.816 


2.461 


1.7-2.322 


0.341 


2.07E+06 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


15 


13.49 


12.54 


0.90 


3.19E-4 


5.293 


2.720 


1.8-2.574 


0.344 


5.66E+06 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


16 


14.16 


11.48 


0.83 


2.78E-4 


5.765 


2.989 


2.0-2.734 


0.347 


6.69E+06 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


17 


14.83 


10.61 


0.76 


2.49E-4 


6.209 


3.253 


2.0-3.151 


0.344 


7.34E+06 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


20 


16.31 


8.68 


0.64 


1.91E-4 


7.643 


4.354 


1.8-4.354 


0.321 


7.08E+06 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


25 


16.35 


6.87 


0.53 


1.40E-4 


10.223 


6.148 


2.2-6.148 


0.287 


6.02E+06 


0.00E+00 


0.00E+00 


0.00E+00 


0.00E+00 


30 


12.92 


5.81 


0.47 


1.01E-4 


12.680 


8.011 


1.9-7.154 


0.263 


5.19E+06 


5.05E+04 


5.05E+04 


0.00E+00 


0.00E+00 


35 


11.94 


5.12 


0.44 


0.80E-4 


14.687 


8.508 


2.3-8.203 


0.251 


4.62E+06 


2.68E+05 


1.16E+05 


1.52E+05 


0.00E+00 


10 


12.52 


4.64 


0.42 


0.78E-4 


16.493 


8.983 


2.5-8.600 


0.247 


4.16E+06 


3.41E+05 


1.10E+05 


2.04E+05 


2.64E+04 


60 


17.08 


3.64 


0.36 


0.62E-4 


25.172 


12.623 


3.7-11.623 


0.218 


3.18E+06 


3.51E+05 


6.18E+04 


1.28E+05 


1.61E+05 


80 


22.62 


3.18 


0.33 


0.57E-4 


34.706 


17.408 


5.0-16.332 


0.187 


2.66E+06 


3.19E+05 


4.43E+04 


8.35E+04 


1.91E+05 


120 


30.83 


2.76 


0.30 


0.50E-4 


48.553 


24.519 


7.8-22.790 


0.156 


2.19E+06 


4.33E+05 


1.78E+05 


5.35E+04 


2.02E+05 


60t 


19.93 


3.93 


0.34 


0.64E-4 


30.034 


15.815 


4.4-14.042 


0.201 


3.20E+06 


3.36E+05 


5.28E+04 


1.02E+05 


1.81E+05 


120t 


26.83 


2.92 


0.31 


0.53E-4 


44.400 


21.360 


6.6-19.601 


0.170 


2.58E+06 


5.35E+05 


2.48E+05 


5.79E+04 


2.28E+05 



OS 



Table 3. 26 A1 and b0 Fe yields 



M ini 


26 Al 


Wind 


Ne/C sh 


cxpl 


60 Fe 


He rsh 


-H-Cg c b:warz. 


-H-CLcdoux 


C s h 


cxpl 


M© 


Mq 


M 


M 


M 


M 


M 


M Q 


M 


M 


M 



11 

12 
13 
14 

15 
16 
17 

20 
25 
30 
35 
10 
60 
80 
120 
60t 
120t 



1.60E-5 
2.11E-5 
2.45E-5 
1.04E-4 
1.32E-4 
6.80E-5 
6.87E-5 
5.43E-5 
8.61E-5 
9.93E-5 
8.38E-5 
1.21E-4 
2.52E-4 
4.00E-4 
7.03E-4 
2.97E-4 
6.97E-4 



1.16E-11 
2.91E-11 
1.46E-10 
4.37E-10 

1.17E-09 
2.74E-09 
6.76E-09 
4.32E-08 
3.93E-07 
2.39E-06 
1.14E-05 
2.06E-05 
6.94E-05 
1.32E-04 
2.82E-04 
5.98E-05 
3.50E-04 



2.20E-6 
2.30E-6 
2.50E-6 
7.00E-5 
8.00E-5 
1.00E-5 
1.50E-5 
1.80E-5 
3.46E-5 
2.61E-6 
2.06E-5 
3.44E-5 
5.06E-5 
8.80E-5 
1.58E-4 
5.02E-5 
1.10E-4 



1.38E-5 
1.71E-5 
2.20E-5 
3.40E-5 
5.20E-5 
5.80E-5 
5.37E-5 
3.63E-5 
5.11E-5 
9.43E-5 
5.18E-5 
6.60E-5 
1.32E-4 
1.80E-4 
2.63E-4 
1.87E-4 
2.37E-4 



1.71E-6 
4.33E-6 
7.56E-5 
5.72E-6 
6.28E-6 
4.39E-6 
7.96E-6 
1.56E-5 
3.69E-5 
1.49E-5 
4.03E-5 
6.23E-5 
2.27E-4 
7.55E-4 
9.93E-4 
5.50E-4 
8.60E-4 



3.0E-7 
2.0E-6 
2.0E-6 
1.8E-6 
8.0E-7 
6.0E-7 
2.0E-7 



— 




1.4E-6 


1.0E-8 







2.3E-6 


7.4E-5 


— 




3.9E-6 





— 




5.4E-6 


8.0E-8 


— 




6.0E-7 


3.2E-6 


— 




4.0E-6 


3.8E-6 


— 




1.4E-5 


1.6E-6 


— 




3.2E-5 


4.9E-6 


— 




7.0E-6 


7.9E-6 


— 




3.3E-5 


7.3E-6 


l.E-5 





4.4E-5 


8.3E-6 


l.E-4 


6.0E-5 


8.0E-5 


4.7E-5 


6.E-4 


1.4E-4 


1.0E-4 


5.5E-5 


8.E-4 


2.5E-4 


1.3E-4 


6.3E-5 


4.E-4 





1.0E-4 


5.0E-5 


7.E-4 





1.0E-4 


6.0E-5 



to 



-63- 



Table 4. Tests on the Al yields 



test 


component 


25 


60 


120 






M 


M 


M 


25 Mg(initial)x2 


wind 




1.34 (-4) 


5.64(-4) 


25 Mg(p, 7 ) 26 AlxlO 


wind 




4.16(-5) 


2.20(-4) 


25 Mg(p, 7 ) 26 Alx3 


wind 




4.60(-5) 


2.33(-4) 


25 Mg(p, 7 ) 26 Alx0.33 


wind 




1.09(-4) 


3.32(-4) 


25 Mg(p, 7 ) 26 Alx0.1 


wind 




7. 79 (-5) 


2.10(-4) 


0.5 H p over, in H burn 


wind 




5.98(-5) 


3.50(-4) 


0.5 H p over, in Hbum. 


total 




2.97(-4) 


6.97(-4) 


26 Al(n,p) 26 Mgx2 & 










26 Al(n,a) 23 Nax2 


expl. 


3.13(-5) 


7.80(-5) 




24 Mg(n, 7 ) 25 Mgx2 


expl. 


8.00(-5) 


2.11(-4) 




25 Mg(p, 7 ) 26 Alx2 


expl. 


8.20(-5) 


2.09(-4) 





Table 5. Selected properties of the models computed with the Langer (1989) mass loss rate 



M ini 


Mend 


M co 


26 A 1 

A1 Wind 


26 A 1 

n - 1 Cshcll+cxpl. 


26 A 1 

A1 Total 


60 w„ 

re Total 


M 


M 


M 


M 


M 


M 


M 



40 


6.88 


5.25 


2.17(-5) 


1.053(-4) 


1.27(-4) 


4.54(-5) 


60 


6.01 


4.54 


7.51(-5) 


4.490(-5) 


1.20(-4) 


1.52(-5) 


80 


6.10 


4.61 


1.42(-4) 


1.120(-4) 


2.54(-4) 


2.17(-5) 


120 


6.18 


4.68 


3.02(-4) 


4.200(-5) 


3.44(-4) 


1.75 (-5) 



